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0. Introduction 

The purpose of this article is to provide an effective method for calculating the wave invariants associ- 
ated to a non-degenerate elliptic closed geodesic 7 of a compact Riemannian manifold (M,g): that is, the 
coefficients in the singularity expansion 

TrU(t) = e (<) + ^e 7 (i) 
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of the trace of the wave group U(t) := expit\fK at t = L 7 (the length of 7.) We will show that 

a l k — / ^ 7 ;fc(s)ds 
J 7 

for certain homogeneous invariant densities / 7 k(s)ds on 7, given by at most 2k+l integrals over 7 of polyno- 
mials in the curvature, Jacobi fields, length, inverse length and Floquet invariants j3j :— (1 — e IQj along 7. 
These expressions characterize the wave invariants in much the same way that the heat invariants (or wave 
invariants at t=0) are characterized as integrals J M Pj(R,V R, ...)dvol of homogeneous curvature polynomi- 
als over M [ABP] [Gi]. Moreover, in combination with the recent inverse results of Guillemin [G.1,2], the 
method produces a list of new spectral invariants of this kind, simpler than the wave invariants themselves 
(the so-called quantum Birkhoff normal form coefficients B^k-.j-) 

To state the results, we will need some notation. We let <g) C denote the space of complex normal 
Jacobi fields along 7, a symplectic vector space of (complex) dimension 2n (n=dim M-I) with respect to the 
Wronskian 

. { X,Y)= 9{ X^Y)- 9 ^X,n 

The linear Poincare map P 1 is then the linear symplectic map on J^-®€- defined by P 7 T(t) = T(t + L 7 ). We 
will assume 7 to be non-degenerate elliptic, i.e. that the eigenvalues of P 1 are of the form {e ±laj ,j — I , n} 
with (Floquet) exponents {ai,...,a„}, together with ir, independent over Q. The associated normalized 
eigenvectors will be denoted {Yj, Yj.j — 1, n}, 

P~/Y 3 = Y J P^Y, = e-^Yj ou(Y 3 ,Y k ) = 5 jk 

and relative to a fixed parallel normal frame e(s) :— (ei(s), e„(s)) along 7 they will be written in the 
form Yj(s) = X)fc=i Ujk{s)ek(s). The metric coefficients will always be taken relative to Fermi normal 
coordinates (s,y) along 7. The mth jet of g along 7 will be denoted by j™g, the curvature tensor by R 
and its covariant derivatives by V m R. The vector fields Jj, and their real linear combinations will be 
referred to as Fermi normal vector fields along 7 and contractions of tensor products of the V m i?'s with these 
vector fields will be referred to as Fermi curvature polynomials. Such polynomials will be called invariant 
if they are invariant under the action of 0(n) in the normal spaces. Invariant contractions against and 
against the Jacobi eigenfields Yj , Yj , with coefficients given by invariant polynomials in the components yjk , 
will be called Fermi-Jacobi polynomials. We will alse use this term for functions on 7 given by repeated 
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indefinite integrals over 7 of such FJ polynomials. Finally, FJ polynomials whose coefficients are given by 
polynomials in the Floquet invariants f3j = (1 — e taj )^ 1 will be called Fermi- Jacobi-Floquet polynomials. 
We give weights to the variables gij, D% y gij, L := L^^ctj^yij^ijij as follows: wgt(D^ y gij) — — \(3\, wgt(L) = 
1, wgt(ctj) = 0, wgt(yij) = |, wgt(yij) = —\. As will be seen, these weights reflect the scaling of these objects 
under g — > e 2 g. A polynomial in this data is homogeneous of weight s if all its monomials have weight s 
under this scaling. 

Theorem A Let 7 be an elliptic closed geodesic with {a\, a n , it} independent over Q. Then a 7 k = 
! 1 I 1 k{s\g)ds where : 

(i) I-yk( s , 9) is a homogeneous Fermi- Jacobi-Floquet polynomial of weight -k - 1 in the data {L, y^ , y^ , D^ y g} 
with |/3| < 2k + 4 ; 

(ii) The degree of L ( k in the Jacobi field components is at most 6k+6; 
(in) At most 2k+l indefinite integrations over 7 occur in J 7 fc/ 

(iv) The degree of 1^ in the Floquet invariants f3j is at most k+2. 

For instance, in dimension 2 the residual wave invariant a 7G is given by: 

c 3 
a 7° ~ ~JJ# [-B 7 o;4(2/3 — (3 — -) + £? 7 o ; o] 

where: 

(a) c 7 is the principal wave invariant i a L#\I — -P 7 | _5 ; 

(b) L # the primitive length of 7; a is its Morse index; P 7 is its Poincare map; 

(c) -B 70 ;j have the form: 

1 f L * 
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0<m,n<3:m+n=3 IV ' Jl J ° 

+ 7^ ]T C 2 , mn Im{f L T„(s)Y m -Y n (s)[[ S T„(t)Y n -Y m (t)dt]ds} 

0<m,n<3;m+n=3 ^ ° J ° 

for various universal (computable) coefficients. Here, 

(d) r denotes the scalar curvature, t v its unit normal derivative, r„„ the Hessian Hess(r){v,v); Y denotes 
the unique normalized Jacobi eigcnficld, Y its time-derivative and Sj Q the Kronecker symbol (1 if j=0 and 
otherwise 0.) 

We note that the residual wave invariant already saturates the description in Theorem A. 

This characterization of the wave invariants makes more concrete (for Laplacians) the recent results of 
Guillcmin [Gl,2] which show that the wave invariants may be expressed in terms of the quantum Birkhoff 
normal form coefficients of the wave operator around 7. For instance, Guillcmin shows [G.2, (8.24)]: 

d 2 11 
a 10 = *c 7 a[]T Qj-Qf H ^ a ^ + 2 )(/3j + 2 } 

+ \ E( 2 & 2 + & - l)|Ui(0, a) - + \)^fM^ °) + 

i 1 i 1 

where i/i_ r (<7, h, . . . , I n ) is the term of order 1-r in the complete symbol of the quantum normal form. The 
coefficients B JO - : j above are essentially these QBNF (quantum Birkhoff normal form) coefficients. 

Indeed, the first step in the proof of Theorem A (which we call Theorem B) consists in explicitly 
constructing the QBNF for >/~K near 7. The method is different from that in [G.1,2] and effectively calculates 
the QBNF coefficients as integrals of Fermi- Jacobi polynomials. It is based to some extent on the construction 
of a complete set of quasi-modes associated to 7 as presented in Babich-Buldyrcv [B.B] although the emphasis 
is on intertwining operators rather than on quasi-modes per se. In the construction of the intertwining 
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operators it also employs several ideas in Sjostrand [Sj], although it does not begin by putting the principal 
symbol in Birkhoff normal form as is done in [Sj] and [G.2]. 

The second step (which we call Theorem C) consists in calculating the wave invariants from the normal 
form. This is possible because the wave invariants are non-commutative residues of the wave operator and 
its t derivatives and hence are invariant under conjugation by unitary Fourier Integral operators (cf. [G.1,2], 
[Z.l]). The residues of the normal form wave group will be easily seen to be polynomials in the QBNF 
coefficients and in the /?/s. Some dimensional analysis of the wave coefficients then leads to the description 
in the statement of Theorem A. 

Guillcmin [loc.cit] has also proved the remarkable inverse result that, conversely, the QBNF coefficients 
can be determined from the wave invariants associated to 7 and its iterates, and therefore are themselves 
spectral invariants. In view of Theorems A-B this gives a list of new spectral invariants, for instance the 
QBNF coefficients B Jo; j above, in the form of geodesic integrals of FJ polynomials, 

Let us now describe the ingredients of the proofs more precisely. Henceforth we will reserve the notation 
7 for a primitive closed geodesic and will denote its iterates by 7™. To introduce the quantum Birkhoff normal 
form at 7 and its role in the calculations, we note first that the wave invariants associated to 7 are determined 
by the microlocalization of A to a conic neighborhood 

(0.1) \y\<e M <£ 

a 

of the cone R + j thru 7 in T*M — 0. Here, (a, 77) denote the symplectically dual coordinates to the Fermi 
normal coordinates (s,y) above. The microlocalization of A to (0.1) is then given by 

A^ := ip(s, D s , y, D y )*Aip(s, D s ,y, D y ) 

where D Xj :— -J^— and where ip(s,a,y,rj) is supported in (0.1) and identically one in some smaller conic 
ncighborhoood . Qften we omit explicit mention of the microlocal cut-off ip in calculations which are valid 
on its microsupport. Under the exponential map 

exp : N 7 — > M 

along the normal bundle to 7, the localization of g, resp. A, to the tubular neighborhood \y\ < e pulls 
back to a locally well-defined metric, respectively Laplacian, on a similar neighborhood in N~. Hence exp 
conjugates A^, to an isometric microlocalized Laplacian on iV 7 , which we continue to denote by A^,. We are 
thus reduced to calculating the wave invariants of a Laplacian on the model space S\ x 1" at the closed 
geodesic 7 = S\ x {0}, where 5^ := R/LZ. For the sake of simplicity we assume the normal bundle is 
orientable but note that the reduction is valid for immersed as well as embedded closed geodesies. 

We now wish to put A.^ into normal form, which is first of all to conjugate it (modulo a small error) 
into a distinguished maximal abelian algebra A of pseudodifferential operators on the model space xl™. 
Roughly speaking, A is generated by the tangential operator D s := -M^ on Sj^ together with the transverse 
harmonic oscillators 

(0-2) I j =I j (y,D v ):=±(Dl+y*). 

In the construction of the normal form, a special role will be played by the distinguished element 

(0.3) K:=j(LD s + H a ) 

where 

1 ™ 

(0.4) H °~oYl akh 

Z fe=i 

and where the choice of sign in ±a/c will be specified below. This element comes up naturally as the semi- 
classical parameter in the construction of quasi-modes, although D s is more suitable for analysing the wave 
invariants. Note that both are elliptic elements in the conic neighborhood 

(0-5) \Ij\<ea Ij(v,v)~\(Vi+T$), 
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which will be the image of (0.1) under the conjugation to normal form. 

The classical Birkhoff normal form theorem states roughly the following: near a non- degenerate elliptic 
closed geodesic 7 the Hamiltonian 



71+1 

H{x,0 = \t\ : X !>'<>^ 
\ ij=i 

can be conjugated by a homogeneous local canonical transformation x to the normal form 
(0.6) X *H^a + ^-J2 a ^ + M \" In) +- mod °lo 

where p k is homogeneous of order k+1 in I\, . . . , I n , and where is the space of germs of functions homo- 
geneous of degree 1 which vanish to infinite order along 7. Note that all the terms in (0.6) are homogenous 
of degree 1 in (a, Ii, I n ), and that the order of vanishing at |/| =0 equals one plus the order of decay in 
a. The coefficients of the monomials in the Pj(h, ■ ■ ■ , I n ) are known as the classical Birkhoff normal form 
invariants. (See the Appendix for some further details). 

The quantum Birkhoff normal form is the more or less analogous statement on the operator level. In 
the following the symbol = means that the two sides agree modulo operators whose complete symbols are 
of order 1 and vanish to infinite order on 7. Also, O m ^ r denotes the space of pseudodifferential operators of 
order r whose complete symbols vanish to order m at (y, n) = (0, 0). 

Theorem B There exists a microlocally elliptic Fourier Integral operator W from the conic neighborhood 
(0.1) o/M+7 m T*N 1 - to the conic neighborhood (0.5) ofT*Sl m T*(Sl x R n ) such that: 

v.-w^a^w =^n,h, ...,/„) [11+ — — — + (L1Z) 2 +■■■+ {Llz)k+ i +•••] 

= n }_tt Pi (h, 1 P2{h,---,In) Pk+l(h,---,In) 

- s+ L a+ LD S + (LD S ) 2 + "' + (LD s ) k+1 + '" 

where the numerators Pj(h, I n ),Pj(h, I n ) are polynomials of degree j+1 in the variables Ii,...,I n , 
where ip is microlocally supported in (0.5), and where W^ 1 denotes a microlocal inverse to W in (0.5). The 
kth remainder term lies in the space (B^tg0 2 (k+2-j)^ 1 ~^ ■ 

The QBNF coefficients will by definition be the coefficients of the monomials in the classical action 
(Ij-) variables in the complete Weyl symbols of the operators Pj(h, ■ • ■ , I n )- As mentioned above, the proof 
of Theorem B gives an effective method for calculating them as integrals over 7 of FJ polynomials. 

The asymptotic relation in the above expansion may be viewed in either of two ways: First, as mentioned 
in the statement of the Theorem, the kth remainder is a sum of terms in A of orders 1,0,..., — (fc + 1) where 
the complete symbol of the term of order 1 — j must vanish to order 2(k + 2 — j). This characterization 
of the remainder will play the key role in the calculation of the wave invariants, since terms in the normal 
form with low pseudodifferential order or with high vanishing order make no contribution to a given wave 
invariant. On the other hand, it may be viewed as a semi-classical asymptotic relation with 1Z playing the 
role of semi-classical parameter; thus the theorem gives a semi-classical expansion for V in terms of TZ. This 
point of view comes up naturally in the theory of quasi-modes associated to 7 : Indeed, consider the joint 
A- eigenf unctions 

(0-7) <P° kq (s,y) :=e fc (s)®7 9 (y) 

with ek(s) := e^ lks and with j q qth normalized Hermite function (q e N"). The corresponding eigenvalues 
of V then have the semi-classical expansions 

/no, X __ 1 gl(g) , P2(g) , 

(0.5) \ kq = r kq H 1 g r ... 

rk 1 r kq 
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where 

1 " 1 

(0.9) r fc , = -(27rA + + -)a,-) 

i=i 

are the eigenvalues of 1Z. Here the index q is held fixed as k — > oo. We recognize in (0.8) the familiar form 
of the quasi-eigenvalues associated to 7 (cf. [B.B., ch.9]); hence the intertwining operator W is the operator 
taking the eigenfunctions (0.7) to quasi-modes of infinite order for A at 7. 

As will be seen in (§4), Theorem B implies that the wave invariants of \/A^ are the same as the wave 
invariants of V. The second main step in the calculation of the wave invariants is then the use of the non- 
commutative residue to connect the terms in the normal form expansion with the terms in the singularity 
expansion for TV exp itV. The main point here is that 

(0.10) a 7fc = res D k e lt ^ := Res s=0 TrD k e lt ^ \fK^,~ 6 ', 

with res invariant under conjugation by (microfocal) unitary operators, and depending on only a finite jet 
of the Laplacian near 7. Hence it may be calculated by conjugating to the normal form, and indeed will only 
depend on a finite part of the normal form. Applying D k and formally exponentiating the terms of order 
< —1 in D s we get 

(0.11) res^(A„/i,...,/n)A fe e^=£=r^ 1 ^ 



LD S 

which suggests that the wave coefficient a^ 7 is the regularized trace of the coefficient of -D s _1 in (0.11). This 
is not clear, even formally, since many of the terms of negative order in D s in the exponent have overall order 
1 as pseudodiffcrcntial operators; but it will prove to be the case. Since e lLDs = I on R/LZ the Fourier 
Integral factor in (0.11) is just e tHa . Regarding the regularized traces of the coefficients, we note that 

(0.12) Tre lH " = e^=i {qk + i )ak 

is well-defined as the tempered distribution 

e i(a k +i0) 

(0-13) T(a) = n£_! .. —r- 

on R^. Since its singular support is the union of the hypcrplancs Zfe m := {(011, . . . , a n ) G R™ : = 2i:m\,T 
has smooth localization to neighborhoods of Floquet exponents (a\, . . . , a n ) which are independent of ir over 
Q. Hence the regularized trace is simply the evaluation of the distribution trace at a regular point. Similarly, 
the coefficient of D s ~ k ~ 1 in (0.12) has a distribution trace of the form 

(0-14) Yl ^k,-i(Qi + l,-,<In + l)e i ^ {gk+i)ak 

= J r k -i(D ai ,...,D a jT(a) 

for a certain polynomial Tk-\- Hence this trace is also a locally smooth function in neighborhoods of 
non-resonant exponents. 

Theorem C The wave invariants are given by 

a kl = J 7 k,-i(D ai ,...,D an )T(a). 



The coefficients of the polynomials Tk-\ arc evidently polynomials in the QBNF coefficients and the 
differentiation process produces polynomials in the /3/s. Combined with Theorem B and a dimensional 
analysis, this proves Theorem A. 

The proofs of Theorems A-C also lead to a somewhat simpler proof of Guillcmin's inverse theorem that 
the classical (in fact the full quantum) normal form is determined by the wave trace invariants for all the 
iterates of 7. We only sketch the proof here, assuming the reader's familiarity with the original proof of 
Guillcmin in [G.2]. The key point is to focus on the Floquet invariants [3j := (1 — e ^)^ 1 for all the iterates 
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7™ of 7, that is the residues (0.10) for t = L, 2L, 3L, .... It follows from the calculations in Theorems A-C 
that the wave invariants are polynomials in the /3j's: more precisely, for each m, a^m is the special value at 
6j = mctj of the fixed polynomial 7 7; fe in the variables /3j. Under the irrationality condition above, the points 
(e mQl , . . . , e miQ ") form a dense set on the torus, and hence the special values at these points determine the 
entire polynomial. The coefficients of the kth polynomial J 7; fc are therefore determined by the wave invariants 
for 7,7 2 , . . . . By studying the relation of the coefficients of I Tt k to the normal form invariants, Guillcmin 
proves that all of the latter can be determined from the former. 

Although Theorems A-C are only proved here under the hypothesis that 7 is non-degenerate elliptic, 
they have analogues for hyperbolic and mixed hyperbolic-elliptic geodesies, which we plan to describe in a 
future article [Z.3]. We also note that for closed geodesies possessing neighborhoods in which the metric 
has no pairs of conjugate points for t < L, the wave invariants can be calculated directly from a Hadamard 
parametrix [D] [Z.1,5]. Since a sufficiently large number of iterates of an elliptic closed geodesies will always 
contain pairs of conjugate points , but a small iterate may contain none, the calculation here and in [Z.l] 
overlap but are independent. The calculations of [Z.l] also apply to hyperbolic geodesies without pairs of 
conjugate points, showing that the form of the wave invariants is essentially the same for the hyperbolic 
and elliptic cases. However the form resulting from the Hadamard parametrix is not immediately that of 
FJF polynomials, and it takes considerable manipulation to show that the formulae given here and in that 
paper agree. In the opposite extreme of Zoll manifolds, all of whose geodesies are closed and of completely 
degenerate elliptic type, the wave invariants are calculated in [Z3,4] by yet another method. 

The organization of this paper is as follows: 

§1: The models 

§2: Semi-classical normal form of the Laplacian 

§3: Normal form: Proof of Theorem B 

§4: Residues and wave invariants: Proof of Theorem C 

§5: Local formulae for the residues: Proof of Theorem A 

§6: Quantum Birkhoff normal form coefficients 

§7: Explicit formulae in dimension 2 

§8: Appendix: The classical Birkhoff normal form 

§9: Index of Notation 

The author wishes to thank S.Graffi for discussions of quantum Birkhoff normal forms during a visit 
to the University of Bologna in June, 1994, where this work was begun. He also wishes to thank Y.Colin de 
Verdiere for his remarks on a preliminary version of the results which were presented during a visit to the 
Institut Fourier in January 1995. The final version was completed after the author received a copy of the 
article [G.2] of Guillcmin, and has benefited a good deal from the discussion there of normal forms. 

1. The models 

As mentioned above, the calculation of the wave invariants associated to a closed geodesic 7 of a 
Ricmannian manifold (M, g) can be transplanted to the normal bundle N 1 by means of the exponential 
map. Thus, the model space is the cylinder S\ x K™, where as above S\ = R/LZ. In this section we first 
collect together some basic formulae and facts concerning the "Hcrmite package" on this model space: that 
is, those aspects of analysis which come from the representation theory of the Heisenberg and metaplectic 
algebras. We then transfer the Hermite package to iV 7 in a way particularly well-adapted to the metric along 
7- 

§1.1: The model: H = H 2 (Sl) <g> L 2 (« n ) 

Since we are only concerned with the conic neighborhood (0.2) of K + 7, we only consider the positive 
part T*Sl x T*R n and its quantum analogue the Hardy space H 2 (Sl) <g> L 2 {« n ) with k > 0. 

On the phase space level, the model is roughly T*(Sj j x R"). More precisely it is the cone (0.1) in 
the natural symplectic coordinates (s,a,y,rj). Since (0.1) is a conic neighborhood of R + 7, we will view it 
as a subcone of the positive part T+S 1 ^ x T*R" (a > 0). On the Hilbert space level the model is then 
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TL := H 2 (Sj^) <g> L 2 (R"), where H 2 {S\) is the Hardy space, or more precisely the range of the microlocal 
cutoff ip of the introduction; generally we omit the subscript unless we need to emphasize the role of tjj. We 
now introduce some distinguished algebras of operators on the model space. 

First is the (complexified) Heiseberg algebra h„®C, which will be identified with its usual realization on 
L 2 (R"). It is then generated by the elements yj = "multiplication by yf and by D yj = j§ y - 1 or equivalently 
by the annihilation, resp. creation, operators 

Aj :=y-j+ iD yj A* = y 3 - iD yj 

which satisfy the commutation relations 

[A j ,A k ] = [A*,A* k ] = Q [A j ,A%]=25 ij I. 

The enveloping algebra of the Heisenberg algebra 

(1.1.1) £:=<Y 1 ,...,Y n ,D 1 ,...,D n > 

is the algebra of partial differential operators on R n with polynomial coefficients. We let £ n denote the 
subspace of polynomials of degree n in the variables yj,D y .. In the usual isotropic Weyl algebra W* of 
pseudo-differential operators on R™, the operators yj,D Vj are given the order |, so that 

(1.1.16) £ n C W" /2 

[£ m ,£ n ] C £ m +"- 2 . 

The symplectic algebra sp(n, C) is represented in £ 2 by homogeneous quadratic polynomials in Yj,Dj, and 
a maximal abelian subalgebra of it is spanned by the harmonic oscillators (0.2). We denote by 

(1.1.2a) I :=< h,...,I n > 

the (maximal abelian) subalgebra they generate in W, with l k :=Ifl W , and by 

(1.1.26) V X = 1C\£ 

the subalgebra of polynomials in the generators (0.2)), with P£ the space of polynomials of degree k. 
The full pseudo-differential algebra on S\ x K™ is the doubly filtered algebra 

x R") = ^*(Sl)(g>W* 

with 

^ mn (sl x r") = o w n . 

A maximal abelian subalgebra of it is given by 

(1.1.3) A:=< D s ,I u ...,I n >=<Tl,h,...,I n > 

where 1Z is the distinguished element (0.3). It inherits a double filtration A mn . As above, our interest is 
really in the microlocalization of (1.1.3) to the cone (0.2), i.e. the operators in (1.1.3) will only be used 
in composition with the microlocal cut-off ip(TZ,I\, ...,/„). In this cone D s and 1Z are elliptic. Hence the 
subalgebra 

(1.1.4) <TZ>(g>Vx 

of pseudo-differential symbols in 1Z with coefficients in V% is well defined. 

An orthonormal basis of L 2 (W l ) of joint eigenfunctions of X is provided by the Hcrmite functions 
j q ,q € N™. Here, -f is the Gaussian 7 Q (j/) = 7i/(y) := e i Ifl . It is the unique "vacuum state", i.e. 
the state annihilated by the annihilation operators. The qth Hermite function is then given by 7, := 
C q A\ qi ...A* n ^ lo {q e N"}), with C q = (2 7 r)-™/ 2 (q!)- 1 / 2 ,a! = qi \...q n \. The notation " 79 " for "Gaussian" is 
standard in this context, see [F], and should not be confused with the notation for closed geodesies. An 
orthonormal basis of _ff 2 (S'| / ) (g> L 2 (R") of joint A- eigenfunctions is then furnished by 



(1.1.5) <f>kq( s >V) '■= e fe( s ) ®7g(2/), e k {s):=e-? 



■iks 
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§1.2: The twisted model H a 

We now introduce a unitarily equivalent (twisted) version of the model, in which the distinguished 
element 1Z gets conjugated to D s . This will eventually help to simplify the transport equations in §2. 
The unitary equivalence will be given by conjugation with the unitary operator 

r© /■© 

'ds 

where \i is the metaplectic representation, and where r a (s) is the block diagonal orthogonal transformation 
on R 2n with blocks 



/■© /■© 
(1.2.1a) Kr a ):= n(r a (s))ds = e^ H " 

Jsl Jsl 



(1.2.16) r aj (s):= 
The direct integral here refers to the representation of i 2 (5' 1 ) ® L 2 (R n ) as J gl L 2 (R n )ds, that is, 

H(r a (s))dsf(s,y) = fi(r a (s))f(s,y) 



f 

Jsl 



where the right side is the application the operator in the y-variables with s fixed. As will be seen below, 
/«(r a ) conjugates 1Z to D Sl and commutes with Ii, . . . , /„. Hence it preserves the algebra (1.1.3). 

On the other hand, it does not preserve the Hilbert space H. Indeed, the elements 4>° kq get transformed 
into the elements 

(1.2.2) e kq (s)(g> lq (y), e kq (s) := e lr ^ s 

which are not periodic in s. Rather they satisfy 

(1.2.3a) e kq <g> 7 (s + L,y) = e tK "e kq ® 7 (s, y) 

with 



n 1 

(1.2.36) N=E(* + 9) a ^ 

The space H 2 (S\) ® L 2 (R") thus gets taken to the space H a of elements of the form 



(1-2.4) f(s, y) = E A*' ?)e fe(Z ® 7 9 

with square summable coefficients. 

We can better describe this Hilbert space (and its associated phase space) in the language of 'quantized 
mapping cylinders'. 

On the phase space level we have the symplectic map r a {L) of T*R n of (1.2.1b), essentially the Poincare 
map of our problem. As in [F.G][G.2] we can introduce its (homogeneous) symplectic mapping cylinder 
Cr a (L) : namely, the quotient of T*R x T*R n under the cylic group < R a {L) > generated by the symplectic 
map 

(1.2.5) f a (L) : T*R x T*R n -» T*R x T*R r \ f a (L)(s,a,y,r]) := {s + L,a,r a (L)(y,r))). 

Note that the first return time is constant, which is consistent with [F.G][G.2, 2.10] since elements of 
Sp(2n, R) preserve the contact form ds— ^(r]dy — ydr]). Also note that the mapping cylinder can be untwisted 
via the symplectic map 

(1.2.7) R a :T*RxT*R n ^T*RxT*R n , R a (s,a,y, V ) := (s, o + \ ^ a M + r« (*)(», v))- 



2 
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Indeed, we have 

(1.2.8) R a (s + L,a,y,r]) = r a (L)R a (s,a,y,r]) 

so that R a induces a symplectic equivalence C r ij\ ~ T*{S\) x R*(R n ). 



On the quantum level, the analogue of the mapping cylinder is the Hilbcrt space H a of functions on 
R x R n satisfying 

(1.2.6) f(s + L,y) = n(r a (L))f(s,y) 

and square integrable on [0, L) x R n . The intertwining operator fi(r a ) is essentially the analogue of R a . 
More precisely, we have: 

(1.2.9) Proposition 

(i) fi(r a )*D s fi(r a ) = TZ; 

(ii) n(r a )<j>° kq (s,y) = e vrk « s lq \ 
(Hi) n(r a ) ■ H — ► H a - 

Proof: 

(i) Follows from the fact that \i takes the jth diagonal block 

1 

to Ij and hence that fi(r a (s)) = e l i Ha . 

(ii) Follows from (i) and the fact that -f q is an eigenfunction of eigenvalue qj + \ of Ij . 

(iii) Follows from (ii). □ 

Under conjugation by /u(r a ), the sub-algebra (1.1.4) goes over to the algebra 

(1.2.10) <D S > 



of pseudodifferential symbols in D s with coefficients in the Ij 's. As usual, it is understood to be microlocalizcd 
to the cone (0.2). 

We now conjugate by a further unitary equivalence to transfer this Hermite package to iV 7 . It will 
induce a new model better adapted to the geometry of (M, g) near 7; we call it the adapted model. 
§1.3: The adapted model 

To define it, we must discuss the Jacobi equation along 7. Let Y(s) be a vector field along 7, and 
as above let Y(s) — Y^i=i Yi( s ) e i( s ) De its expression in terms of the parallel normal frame. The Jacobi 
equation is then 

d 2 " 

_y,+ £^.y, = 0. 

ij=l 

Let Jj- ® C denote the space of complex orthogonal Jacobi fields along 7. Equipped with the symplectic 
form lo(X,Y) :=< X, > — < ^r-,Y >, it is a symplectic vector space of dimension 2n (over C). Here 
^ denotes covariant differentiation along 7, and <,> is the inner product defined by the metric. Now 
let Z(s) := (Y(s),^f). The linear Poincarc map P 7 is by defintion the operator on ® C given by 
PjZ(t) = Z(t + L). We recall that 7 is assumed non-degenerate elliptic, hence 

Spec{P 7 ) CS 1 ± 1 i spec(P 7 ). 

Since P 7 e 5 f p(2n,R), its eigenvalues come in complex conjugate pairs {e Mfe ,e _Mfc }. The exponents au 
(called Floquet) are only defined up multiples of 27r; they will be normalized below as in [B.B, (9.3.17)]. 
Then there exists a basis of complex eigenvectors {Yi, Y n } satisfying 

(1.3.1) P- t Y k = e ia "Y k Yj e J® C 

u(Y k ,Yj) = u;(Y 3 ,Y k ) = 0, u{Y 3 ,Y k ) = 6 jk . 

with one choice of eigenvalue from each complex conjugate pair. Equivalently, the span of {Y\, Y n } defines 
a P 7 -invariant positive Lagrangean subspace of ® C. 
Consider now the (modified) Wronskian matrix 

_ / ImY( S y ImY( S y 
[i - 6 - Z) Gs - { ReY(s)* ReY( S y 
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Here, with a little abuse of notation, we denote the components of Y(s), rcsp. relative to the parallel 
normal frame e(s) by Y(s), resp. Y(s) and the notation A* refers to the adjoint of a matrix A. From (1.3.1) 
we have that a s <E Sp(2n, R) , which is equivalent to 

as ds 

As above, we let \x denote the metaplectic representation. Identifying a s with one of its two possible 
lifts from Sp(2n, R) to Mp(2n, R) we introduce the unitary operator 

/■e 

(1.3.3) n(a) := / p{a s )ds 

J -y 

on J®! L 2 (W l )ds. In other words, 

fJ-(a)f(s,y) = n{a s )f(s,y) 

where the operator on the right side acts in the y-variables. Informally, we think of the range as the Hilbert 
space J L 2 (A^ 7 ( s ))c?s, and below we will describe it more completely in terms of quantum mapping cylinders. 
We now use n(a) to transfer the Hermite package to 7V 7 . We begin with the generators of the Heisenberg 
algebra, and set: 

Pj := /z(a)* Dj/i(a) Qj := /j,(a)*Yjjj,(a) 
Aj := n{a)*Ajn(a) A* = p(a)*A*n{a) 

We will refer to the Aj's, resp. A*'s, as the adapted annihilation and creation operators and to the 
operators J 7 j in the symplectic algebra as the adapted action operators. These adapted operators play a 
key role in the study of quasi-modes associated to 7. To establish the connection, we now verify that they 
coincide with the operators similarly denoted in [B.B., ch.9]. 

First, some remarks on notation. For reasons that will become clearer below (§1.4-5), we change the 
notation for the transverse coordinates from y to u, which should be thought of as the rescaled coordinates 
u = L~?y. Objects in the adapted model will henceforth always be expressed in terms of u— coordinates. 
For instance, multiplication by Uj will be denoted simply by uj and differentiation in Uj by D Uj . Also, to 
quote easily some basic facts about the metaplectic representation from Folland [F], we will conform to the 
following 'transposed' notation for the remainder of this section: symbols of operators will be denoted a(p, q) 
rather than a(q,p), and the corresponding Weyl pseudodifferential operator will be denoted by er(D, x) (later 
we will also use the more standard notation a w (x,D)). An element of Sp(n, K) will be denoted by A and 
its transpose by ^4* . 

(1.3.4) Proposition 



Proof: This follows from the metaplectic covariance of the Weyl calculus, that is from the identity [F, 
Theorem 2.15] 

{a- A){D, x) = niA^aiD.x)^*)- 1 . 
If we set a = A* we get (in an obvious vector notation) 

li(a)* (uj + iD Uj )fi(a) = i(Yj ■ D u ) - Yy u 
(compare [BB, ch.9]). □ 
As for the tangential operator, we have: 
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(1.3.5) Proposition The image C of D s under [i is given by: 

^ n n 

C := ^a)*D s ^a) = D s - - D 2 U , , + £ K t] { S )u lUj ). 



2 V 

Proof: The left side is equal to (D s + fi(a s )* D s [i(a s )). To evaluate the second term, we use that both 
ReY(s) and ImY(s) are Jacobi fields, and that Jacobi's equation is equivalent to the linear system jj^(Y, P) = 
JH{Y, P). Here, P — J is the standard complex structure on R 2n , and 

H 



K 
/ 

where K is the curvature matrix and / is the identity matrix [B.B, (9.2.9)]. Hence, the second term is 
ld/j,(JH) with dfj, the derived metaplectic representation. But id/j,(JH) = 1/2(^™ =1 d 2 . — 5Z™- =1 Kij(s)uiUj) 

h ~ '' □ 

We now consider the appropriate Hilbert space (quantized mapping cylinder) in the adapted model. 
We first note: 

(1.3.6) Proposition 

(i) 

M ( a - 1 ) 7o (s, U ) := U (s,u) = (rfetF( S ))- 1 / 2 e i 5<r( S )« > «> 

where T(s) := %Y^. 
(ii) 

M (a- 1 ) 79 :-C/ 9 = A? 1 ...A^[/ . 



Proof: Let us recall the action of the metaplectic representation on Gaussians [F, Ch. 4.5]. For Z in the 
Siegel upper half space (n x n complex matrices Z = X + iY with Y » 0), we have the Gaussian 

lz (x) : = S <Zx ' x> 

on W 1 . The action of an element A <G Mp(2n,M.) on the Gaussian is given by 

M(^* _1 )7z = m(A,Z)~/ a(A)z 

where 

m(A, Z) = det- 1/2 (CZ + D),a(A)Z = (AZ + B){CZ + D)^ 1 

for 

A=( A B 
\C D 

(see [F, 4.65].) Writing 

Y(s) = ImY(s)i + ReY(s), T(s) = {ImY(s)i + ReY)(ImY(s)i + ReY(s))- 1 

we see that m(A, H)la(A)u{ u ) — {detY{s))^ 1 / 2 exp{^ < T(s)u,u > if 

. = / ImY(s) ReY \ 
\ ImY(s) ReY(s) J ' 

The formula (i) follows from A = a* and (ii) is an immediate consequence of (i). □ 
Consider now the periodicity properties of the above data under s — > s + L. We observe that a s fails 
to be periodic in s for two reasons: first, due to the holonomy of the frame e(s), and secondly due to the 
monodromy of the Jacobi fields Y(s). Indeed, we have: 

(1.3.7a) a(s+L) = Tasr' 1 ^) 



with r a as in §1.2 and where T, the holonomy matrix, is the 2n by 2n block diagonal matrix with equal 
diagonal blocks t — (Uj) satisfying 

n 
3 = 1 

It is of course the lift to T*R™ of a rotation on the base. The two properties can be summarized by writing 
(1.3.76) 5 S+L =T5 S 
with a s :— a s r a (s). 

As in §1.2 we will reformulate (1.3.7a-b) in terms of quantum mapping cylinders. First, we put 
(1.3.8) C??(R x R") := {/ g C°°(]R x R n ) : f(s + L,u) = /i(T)/(s,u)} 

and let Ht denote its closure with respect to the obvious inner product over [0, L) x R". Note that the 
metaplectic operator n(T) is simply 

MT)/(u) = fit-'u) 

and hence that 

CfP(R x R n ) ~ C°°(iV 7 ) 

where the isomorphism is simply the pull-back by the exponential map defined by the frame e(s). Thus: 

(1.3.8) Proposition 

(i) Let P(s,D s ,u,D u ) be a partial differential operator on Nj, expressed in the coordinates (s,u). Then 

P(s + L, D s ,u, D u ) = n{T)P{ 8 , D s ,u, D u )n{T)*; 

(ii) The functions 

/i(5 a )($,)=&, :=e^ s U q (s,u) 
define a smooth orthonormal basis ofTir- 
Proof 

(i) It suffices to prove this when P is a vector field given in the local normal coordinates by a (s, u)D s + 
Y^j=\ o-j{s,u)D Uj . Since the metaplectic operator n{T) corresponding to T is the operator f{u) — > /(t _1 ii), 
we have 

n 

n(T)Pn(T))* = a (s, t~ 1 u)D s + ^ a {s, t~ 1 u)t jk D Uk . 

3 = 1 

On the other hand, the vector field is well-defined on iV 7 if and only if 

a Q (s + L, t~ 1 u)D s + ^2 a i( s + L, t~ 1 u)ti D U] = a (s,u)D s + ^ aj(s,u)D Uj . 

(ii) Clear, since by (1.3.7b) /x(a) intertwines the model and the quantum mapping cylinder of /i(T). □ 
Remark: Statement (ii) is equivalent to 

Uqis + L^u) = e- lK "U q ( Sl u) 
(correcting the formula stated in [B.B., (9.3.25)].) 

It follows that we may write a smooth function / in the adapted model in the form 

oo 

(1.3.9) f(s,u) =J2H f(k,qy^ s U q (s,u). 

k=o <jSN" 

§1.4 Metric scaling and weights 

As mentioned above, A and hence the wave invariants have well-defined weights under the metric 
rescaling g — > e 2 g. Since the wave invariants will be expressed in terms of QBNF coefficients, it is natural to 
ask how the latter scale. The question is not really well-posed since the QBNF coefficients are coefficients 
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with repsect to Harmonic oscillators d% + y? whose scaling behaviour depends on the choice of coordinates. 
To amplify this point, we record how various metric objects scale under metric re-scaling. 

In the following table, (s, y) denote the Fermi normal coordinates relative to g, (s, u) denote the scaled 
coordinates (s, L~?y) and p u denotes the symplectic coordinates dual to u. 



9 


e 2 9 


L, (s,y) 


eL, (es,ey) 






Vij :=5(^i,ej) 


&vn = e 2 g(€-^Y j ,e- 1 e j ) 


Kij = g(R(d s ,ei)d s ,ej) 




u = L~2-y Pu ^L^rj 


6 2 u e~zp u 


a j = 'E'k=i( i yjkD Uk - yjkUk) 


Y^'k=i(WjkD Uk - y jk u k ) 


£:=D S - ±{YU D l + E«=i K t] { S )u t u 3 ) 





The entry j/y — > t 2 yij follows from the scale invariance of the Jacobi equation together with the 
normalization condition 

g(ReYi, ImYi) — g{ReY il IvnYi) = Constant. 

which implies that Yj — ► e~^Yi. 

We observe that the creation/annihilation operators, hence the harmonic oscillators A* A.,-, of the 
adapted model are scale- invariant, and that the distinguished element C has weight -1. These are the 
desired scaling properties and we would like the basic and twisted models to possess them as well. As they 
stand, these models do not scale properly if we interpret the (s, y)-coordinates as Fermi normal coordinates. 
However, they do scale properly if we interpret the y coordinates as weightless. 

§1.5 Scaled adapted model and intertwining operators 

To avoid confusion, we now introduce the weightless coordinates x — L~ x y = L~?u, £ = Ln = L^p u and 
henceforth use them exclusively for the scaled adapted, basic and twisted models. In the following table we 
record how the various objects appear in the weightless coordinates. We also record the various intertwining 
operators, since they get altered when we used weightless coordinates. For instance, intertwining by fi{r a ) 
above is weightless but that by fi(a) is of mixed weight. 



Adapted model 


Scaled adapted model 




T,k=i(iyjkL-?D Xk -L-i^x k ) 


C 


D s - £(£;=i L-'Dl + £" ?=1 LK ijUiUj ) 


U (s,u) 


U Lo (s,x) = (dety(s))-V2 e -r<r(^>|dx|5 



The following intertwining operators will arise in the construction of the normal form: 



From * to * 


Classical 


Quantum 


Ad. model to Sc. ad. model 


(u,p u ) —> (L-?u,L?p u ) = (x,£) 




Sc. Ad. Model to Tw. model 


-» (L-*ReYx + L-ilmYZ,LiReYx + L*ImY£) 


U ( D L i ■ a ) 



Above, the notation u(D r ) refers to the metaplectic (dilation) operator corresponding to the symplectic 
matrix 

D r := ( q ) , li(D r )f(s,y) := f{8,r^y). 
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2. Semi-classical normal form of the Laplacian 

We return now to \/A, which as in the introduction will be identified with its transfer to N 7 under the 
exponential map. The finite jets of this transfer are globally well- defined on 7V 7 , so we will often treat A as 
if it too were globally well-defined. Our purpose is to define the semi- classically rescaled Laplacian Ah and 
to put A h into a semi- classical normal form. This is the crucial preliminary step in putting A itself into 
normal form. 

In view of §1.5, there are two rescalings at hand: the semi-classical rescaling Ah and the metric L- 
rescaling above. The two rescalings have quite distinct origins, so we have kept them separate. 

To motivate the rescalings and the emergence of semi-classical asymptotics, let us recall that the quasi- 
modes associated to 7 have the form 

00 

(2-1) * kq (s, jrT q y) = e^^r^uhs,^^) 

1=0 

with U° = U q (sec [B.B]). The intertwining operator W 1 to the normal form is then the operator defined by 
the equations 

(2-2) W 1 <f> kg (s,y) = $ k q(s, y /r k - q ~y). 

The higher order terms hence W 7 , are determined by the conditions 

(2-3) A y e ir ^ s U kq {s,^r k -y,r^) ~ \ kq e ir "« s U kq {s, Jrj^y, r^ 1 ) 

with X kq given by (0.9). 
Now write 

1 ™ 1 

Tkq = h^L 1 hkq := ^ + + 2 aj ^~ 1 

so that the Planck constants h kq are metric- independent. In the scaled Fermi coordinates u of §1.3-5, the 
quasi-modes have then the form 

! j a I 

(2.1') ®kq(s,h *u) = e h "" L U kq (s, y/h kq 2 u,h kq ) 

and the eigenvalue problem (2.3) becomes 

(2.3') A u e^ s U(s, h~}u, h kq ) = \{h kq )e^ s U(s, h' k }u, h kq ). 

We are thus led to study the asymptotic eigenvalue problem 

(2.4) A u e^ s U(s, hT*u, h) = \{h)e^ s U{ Sl hr*u, h) 

on C°°(R 1 x W l ) with U(s,u, h) and X(h) asympotic series in h. Since A u comes from an operator on 7V 7 , 
the eigenvalue problem is taking place on the quantized mapping cylinder TCt (§1-3) of the adapted model. 

We also note that the local expression for the Laplacian in the (s,u) coordinates is the same as in the 
(s,y) (Fermi) coordinates: 

1 " 

A « = -7= E ^9 v Vddu, - W Uj ., =L -i V] 

VS ij=o 

since — > L~ 1 g %: > and d Ui — ► L^d Vi . Momentarily we are going to rescale the coordinates again to the 
weightless rr-coordinates of §1.5, and again the Laplacian will be given by the usual expression. Hence it will 
be a simple matter to pass back and forth between the (s,y), (s,u) and (s,x) expressions. 

It is natural at this point to introduce the unitary operators Th and Mh on Ht or cquivalently on the 
1/2-density version L^R 1 x W l ,il 1 / 2 ) given by 

(2.5a) T^/Mldsl 1 / 2 !^! 1 / 2 ) := h~ n ' 2 f{s, h^u^ds^' 2 ^ 1 ' 2 



(2.56). 



M h (f(s,u)\ds\^ 2 \du\^ 2 ) := e^ s fis^ds^duf' 2 
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We easily see that: 

(2.6) T* h D Uj T h = h-i Duj 

T,*UiT h = h^u t 
M* h D a M h = {{hL)- 1 +D S ) 
[M h , Ui ] = [M h ,D Ui ] = [M h ,T h ] = [T h ,D s ] = 0. 
Definition The reseating of an operator A u = a(s, D s ,u 7 D u ) of the adapted model is the operator 

(2.7) A h := T* h M* h AT h M h 

We observe that the operation of rescaling is weightless, fn particular, the rescaled Laplacian in the 
sense of (2.7) is of weight -2. To calculate it, we first note that the (1/2-density) Laplacian in scaled Fermi 
normal coordinates is given by the expression 

n 

(2.8) A u = -{J- 1 d s Jg°°d s + J'^u^Jd^) 

ij=l 

where d x := and where J = J{s,u) = ^fg is the volume density in these coordinates. To obtain a 
sclf-adjoint operator with respect to the Lesbegue density |rfs||du|, we replace A by the unitarily equivalent 
1/2-density Laplacian 

A 1/2 := J^AJ- 1 / 2 , 

which can be written in the form: 



/2 



(2.9) -A 1/2 = J-^dsg^JdsJ- 1 / 2 + J^du^JdujJ- 1 

ij=l 

n n 

= g°°d 2 s + r°d s + 9 ij d Ui d Uj + ]T + °o- 
ij=l i=l 

From now on, we will only use Ai and denote it simply by A. 
We then have: 

(2.10) -M*AM h = -(hL)- 2 g°° + 2i(hL)- 1 g 00 d s + {{fiL^T + A 
Conjugation with Th then gives 

n n 

(2.11) -A h = -(hL)- 2 g$+2i(hL)- 1 g$d s + i(hL)- 1 r° [h] + 

the subscript [h] indicating to dilate the coefficients of the operator in the form, fh(s,u) :— f(s,hiu). 
Expanding the coefficients in Taylor series at h = 0, we obtain the asymptotic expansion 

CO 

(2.12) A, ~ J2 ^ 2+m/2) A- m /2 

m=0 

where £2 = L , £3/2 = an d where 

f) 1 n n 

(2.13) . & = lL-\i- + -{^X - Y, K ij(*)Wi)] 

We observe that L\ is 2L~ 1 times the distinguished element C of the adapted model, and that C has weight 
— 1 when the u-variables are given their natural weights 1/2. 

15 



As discussed in §1. 4-5, it will be helpful to rescale the variables once again to make them weightless. 
Hence we change variables to x = L~^u and rewrite Ah and the 's in terms of the ^-variables. For 

instance we will henceforth write C in the form: 



o -i n n 



ds 2 1 . 

We further note that the operators £ 2 - m / 2 now satisfy the periodicity condition (1.3.8ii) : Indeed, 
as noted above, A has this property when expressed in normal coordinates, and the various transversal 
rescalings and the conjugations by T h ,M h preserve it. To indicate that an operator has this periodicity 
property and hence acts on Ht, we will subscript the appropriate spaces of operators with a "T." 

From (2.11) we see that the terms in C 2 -m arc of the form 

x m x m ~ 2 x m ~ ^ 

x m ~ 2 D 2 x x m - 3 D x 

x m ~ 2 D s x m - 4 D s x m - A D 2 s 

hence 

(2.14a) £ 2 _ f e^fSV^r' + ^^T^r' + ^TX^ 

where the subscript e indicates that the Weyl symbol is a polynomial with the parity of m. Moreover, in the 
s variable, £2-^ is a polynomial in D s of degree at most two, so we can refine (2.14a) to the statement 

(2.146) £ 2 _~ e C^{Sl£T- A )D 2 s {Sl,£?~ 2 )D a +C?{Sl,£™) 

Comparing with (1.3.4) we see that 

C = v(Al)DMAi)- 1 
where Al is the weightless Wronskian matrix D^ia*, that is 



A L ■= 



LzImY L^ReY 
L-^ImY L'^ReY 



This motivates the conjugation of (2.11) to the (untwisted) model. We therefore put 

V h = n(A* L r 1 A h n(A* L ) 

which has the asymptotic expansion 

00 

(2.15) V h ~ J2 h { - 2+ ^V 2 _^ 

m—o 

with T>2 = /,£>! = 0,T>i = D s . Conjugation by n(A* L ) preserves weights, homogeneity and parity in the 
variables (x, D x ). It also transforms D s into D s plus a term quadratic in (x, D x ). Hence we find easily that 
T> 2 -m. has weight -1 for each m and that 

(2.16a) V 2 -m. e * 2 (R) <g> B™- A + *^R) <g> £ £ m ~ 2 + *°(R) ® £™ 

or, analogously to (2.14b), 

(2.166) £> 2 _m g C°°(R,£™- 4 )D 2 S +C°°(R,£™- 2 )D S + C°°(M,£ e m ) 

Of course, conjugation with fi(A* L ) also alters the periodicity property of the terms T) 2 -m.. From 
(1.3.7a) and (1.3.8) we see in fact that they transform like operators in the twisted model, i.e. on the 
quantum mapping cylinder of r a (L). More precisely, we have 
(2.17a) 

2? 2 _ f \ S+L = n(A* L ) ~ 1 (£2- m v(A* L )\ ( S +L) ) = niraWMAiy^iTyCi-y \ {s+L) n(T)(,(A* L )t,(r a (L)y 
= fJ,{r a (L))n(A* L )*C 2 -^n{A* L )fi{r a (Ly = /i(r„(L))D2-M/f(r a (L))*. 
Equivalently, in terms of the matrix elements in the basis of Hermite functions, we have 
(2.176). (V j7q , lr )(s + L) = e-^-^ s {V jlq , lr ){s) 

16 



To indicate that these operators act on H a we henceforth subscript the appropriate spaces of operators with 
an a. 

To render these terms periodic in s, we have to conjugate further to the (untwisted) model under n{r a ). 
We record the resulting expressions, since they will be used later on. In the notation A* L (s) := A* L (s) ■ r a (s) 
the principal operator becomes, by Proposition (1.3.8), 

(2.18) ^{ADCuiA)- 1 = n 

Since conjugation by fi(r a ) also preserves weight, homogeneity and parity, we further have: 

oo 

(2.19) Tl h := n(a)A h »(a)* ~ ]T h^+^K^ 

m—o 

with 7?-2 = I,TZ-3 = 0,7?.i = 1Z and with all coefficients of weight -2 and periodic, that is, 
(2.20a) ^ 2 _™ G ^(S 1 ) <g> £ e m ~ 4 + ^(S 1 ) ® £™~ 2 + ^"(S 1 ) <g> £ £ m 

or, analogously to (2.14b), 

(2.206). K 2 -m. € C°°(,Si,£ £ m - 4 )ft 2 + C°°{S 1 L ,£™- 2 )1l + C°°(S 1 L ,£™) 

Our aim is now to put A, or for certain h, into a semi-classical normal form. This normal form will 
be, at first, only a formal normal form for A^ as a formal ft,— pseudo-differential operator. Here, a formal 
/i-pseudodifferential operator of order m on is the Weyl quantization 

(2.21a) a w (x,hD-h)u(x) = (2irh)- n J I ' e H*-v<$ a (^(x + y),t;h)u(y)dyd£ 

of an amplitude a belonging to the space S™(T*R N ) of asymptotic sums 

oo 

(2.216) a(x,£,h) ~ h- m ^2aj(x,0h j 

j+o 

with a,j € C°°(R 2N ). Such operators form an algebra ^(R^) under composition, with ^™ the subspace of 
mth order elements. (Sec [Sj] for further background and references). We will also be concerned with the 
slightly different situation of h- pseudodiffcrcntial operators on S 1 x R n , which arc defined similarly using 
local coordinates. Combining the h— filtration with the previous nitrations of **(K) <g> W*(R") we get the 
triply filtered algebra 

#(*.*.*) ( K x R n) 

with fc + m + ^ the total order of an element (and similarly with S* 1 replacing R.) 

The following lemma will prepare for the normal form. We state it in terms of the ^-operators of the 
model since the periodicity properties are simplest there. The notation U A\ " will be used for the restriction 
of an operator A G , 5*(5 1 x R") of the model to elements of 7?.-weight zero. Equivalcntly, after conjugation 
by n(r a ) to the twisted model, to elements of D s -weight zero , that is, to functions independent of s in H. a . 
If we write the latter A in the form A 2 D 2 + A\D S + A Q , then A\ a is A \ (restricted to weight elements). 
We have: 

(2.22)Lemma There exists an L-dependent h-pseudodifferential operator Wh = Wh{s, x, D x ) on L 2 (S^ x R n ) 
such that, for each s G S\, 

W h {s,x,D x ) : i 2 (R") ^ L 2 (R") 

is unitary, and such that 

oo 

W* h TZ h W h ~ -h- 2 L- 2 + 2h- 1 L- 1 TZ + J2 hi K™_j_(s, D s , y, D y ) 

3=0 

where 

(i) (s, D s , x, D x ) = 11°°'* TZ 2 + 1l°°')ll + 11°°'° , with 11°°* e C°°(Sl £r 2k ); 

z 2 Z 2 Z 2 Z 2 z 2" 
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(ii) RgLj^jDsjXjDJlo = 7?.~'°(s,a;,£> x )| = fj{h, ■ ■■,I n )\o for certain polynomials fj of degree j+2 
on E™, i.e. fj(h,...,I n ) eVt 2 

(iii) K™_ 2J ^(s,D s ,x,D x )\ = K™j°^{s,x,D x )\ = 0; 

(iv) The harmonic oscillators are weightless and all of the IZ's have weight -2. 
Proof: The operator Wh will be constructed as the asymptotic product 

(2.23) W h := nirayiL^W^nira) 

of weightless unitary /i-pseudodiffcrential operators on E™, with 
(2.23a) W h t := exp(ih^ Q k) 

and with h? Qk e h* C°°(Sl) <g> £ k+2 of total order 1. The product will converge, for each s, to a unitary 
operator in \I/£(E n ) (see [Sj] for discussion of asymptotic products). 

To see what is involved, we first construct a weig htless Qi (s, x, D x ) e C°°(5i) <g> £\ such that 

(2.24a) e~ ihiQ hTZ h e ihhQ i \ a = [-h^L' 2 + Ih^L^ll + K*+...]\ 

where the dots . . . indicate higher powers in h. The operator Qi then must satisfy the commutation relation 
(2.246) {[L- 1 K,Qi]+Ki}\ o = 0. 

To solve for Qi , it is convenient to conjugate back to the T>2-m.'s of the twisted model by /i(r Q ), since this 
transforms 1Z into D s . The commutation relation thus becomes 

(2.24c) {[L- 1 D s , f i(r a )*Qin(r a )}+Vi}\ o = 0, 

that is, 

(2.24d) £ -1 a {Mra)*QiM(ra)}|o = ~i{^i}\o 

where d s A is the Weyl operator whose complete symbol is the s-derivative of that of A. Since (2.24d) is 
simpler than (2.24b), we henceforth conjugate everything by /u(r a ), and relabel the operators n(r a )*QfJ,(r a ) 
by Q. The resulting D's then have the twisted model periodicity properties (2.17a-b). Our problem is then 
to solve (2.24d) with an operator Qi satisfying (2.17a-b). 

To do so, we first observe that under conjugation by /x(r Q ), elements of 1Z- weight zero transform to 
elements of D s - weight zero, and hence it suffices to solve for the matrix elements (Qij qi j r )- It follows from 
(2.17b) that the solution is unique for r ^ q, while for r = q the matrix element (Qi7 9 ,7g) is periodic and 
hence a necessary and sufficient condition for solvability by an operator with the correct periodicity is that, 
for all q G N™ , 

(2-25). J (V i7q , lq )d S = 

In view of (2.16a-b), we have 

Vi e x f ,1 (M) a O £) + *°(R) a ® £?, 

in fact a simple calculation shows that Vi G V I /0 (R) Q <g> £f . (Recall that the subscript a indicates the 
periodicity property (2.17a-b).) 

We now observe that if A e **(R) <g> £™ then 

Li — > L\ m even 

± T 

where (resp. L 2 _) denotes the subspace of L 2 (R) <g> L 2 (R") spanned by even (resp. odd) functions of 
y e E" with coefficients in s. 
It follows that 

Vi : L\ -» L 2 

2 
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(2.26) A : ( L ± ~" "± 

y ' \ : Li —> Li m odd 



This parity reversing property implies the vanishing of the diagonal matrix elements in the basis {7^}, hence 
(2.25) does hold. A unique solution of (2.24d) is specified by the condition that 

{Qi~1 q ,~1 q )ds = 0. 



/ 

Jo 



To solve the equation (2.24d) we rewrite it in terms of complete Weyl symbols. We will use the notation 
A(s,x,£) for the complete Weyl symbol of the operator A(s, x, D x ). Then (2.24d) becomes 

(2.27a) L^dsQ 1 (s, x, £) = -iVi \ a (s, x, f ) 

with 

Qi(s + L,x,£) = Qi{s,r a (L)(x,£)). 
We solve (2.27a) with the Weyl symbol 



Qi(s,x,£) = Qi(Q,x,£) + L / —iVi\ (u,x,€)dv 
2 2 Jo 2 



'0 

where Qi(0, x, £) is determined by the consistency condition 

(2.276) Qi(L,x,£)-Qi{0,x,£) = L [ -iT>i\ (u,x,£)du 

2 2 Jo 2 

or in view of the periodicity condition in (2.27a), 

(2.27c) Ql(0,r a (x,S))-Qi(0,x,£)=L [ -iDi\ (u,x,Qdu. 

Jo 

To solve, we use that T>i\ {u,x,t;) is a polynomial of degree 3 in (a;, £). Also, as is customary in such 
calculations, we switch to complex coordinates Zj = Xj + i£j and Zj = Xj — i£j in which the action of r a (L) 
is diagonal. With a little abuse of notation, we will continue to denote the Weyl symbols, qua functions of 
the Zj,Zj : s by their previous expressions. We also suppress the subscripts by using vector notation z, z and 
e ia . Thus, (2.27c) becomes 

(2.28) Qi (0, e ia z, e~ ia z) — Qi (0, z,z) = L j -iVi\ (u,z,z)du. 

2 2 Jo 2 

We now use that T>i (u, z, z) is a polynomial of degree 3 to solve (2.27c). If we put 

Q1(S,Z,Z)= Yl %mn(s)z m Z n 
|ro| + |n|<3 

and 

V h \ {s,z,z)du= d hmn (s)z m z n 

|m| + |n|<3 

then (2.28) becomes 

(2.29) . V (l-c(™-") Q )gi ;m „(0)z m z" = -zL 2 V di 



|m| + |n|<3 |m| + |ra|<3 



Since there are no terms with m = n in this (odd-index) equation, and since the aj's are independent of it 
over Z, there is no obstruction to the solution of (2.29). 

For simplicity of notation we will express the solution in the form 

(2.30) Qi\o = -iL J Vi(s)ds\ 

where J denotes the integration procedure just defined, that is, the indefinite integral satisfying (2.27b). 
Since the integration only involves the s-cocfficicnts, the solution is a pscudodiffcrcntial operator on K™ with 
the same order, same order of vanishing, and same parity as the restriction of Di to elements of weight zero. 
We then extend it all of to H. a as a pseudodifferential operator with the same properties by stipulating that 

(2.31a) QiM h = M h Qi. 
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Then, as desired, 

(2.316) Qi e *°(M 1 ) ® £ 3 . 

The conjugate by /x(r a ) then defines a periodic operator satisfying (2.24b) and hence a unitary element 
W h ± £ ^hi^ 1 x ^™) satisfying (2.24a). The twisted unitary operator with exponent Qi, i.e. the image of 
W h i under conjugation by /i(r Q ), will be decorated with a tilde, W h i. 

Since h^Qi is of total order 1, h^ad(Qi) (with ad(A)B := [B, A}) preserves the total order in 
and hence W h i is an order-preserving automorphism. It is moreover independent of D s and has an odd 
polynomial Weyl symbol, so that 

(2.32) h^ad{Qi) : h?¥(R)®£™ -» [tf'- 1 ^) ® £ £ m+3 + *'(R) ® £ £ m+1 ]. 

Finally, since 2?± has weight -2, L has degree 1 and since J^(-)ds has degree 1 in the metric scaling, we see 

that Qi is weightless. Alternatively, the di. mrl 's have weight -2, the variables z have weight and hence the 

<7i;m,™' s have wei g ht 0- 

Consider now the element 

-Dl := W* hi V h W h i G ^(M 1 x R»). 
Using only the ^-filtration, we expand 

CO .j 

(2.33) ^-E^ 2+t E 7Wi) 3 ^-f 

n—o j-\-m—n ' 



h- 2 L- 2 + h^L^D, + h- 2+ %V. 



2- 

n— 3 



An obvious induction using (2.32) and (2.16a-b) gives that 

ad(Qi) j V 2 _ f e C°°(R,£^- A )D 2 S + C°° (M.,£™ +j ~ 2 )D s + C°°(K.,£™ +j ). 
i 

It follows that T>2_n. has the same filtered structure (2.16b) as T>2-%. 

We carry this procedure out one more step before arguing inductively, since the even steps behave 
differently from the odd ones. We thus seek an element Qi(s, x, D x ) e 5'*(5' 1 x K") and an element 
fo(h, —,I n ) e Aso that 

V\ := W* hl V^W hl = h- 2 L- 2 + h~ l L- l D s + h~^V\ + V\{s, D„,x, D x ) + . . . 

2 

with 

(2.34a) T>l(a,D a ,x,D x )\ = f (h, I n ) 

with Whi — e 1 ^ 1 , and where the dots signify terms of higher order in h. Note that T>\ = T>1, so that 

2 2 

(2.33a) would imply 

(2.346) {hiv\+Vl}\ = / (Ji, ..,/„). 

2 

The condition on Q\ is then 

(2.35a) {[D s ,Qi]+vl}\ = f (h, I n ) 

or equivalently 

(2.356). d s QA = {-VI + f (h, J n )}| 

As above, we first consider the solvability of (2.35b) one matrix element at a time. The off-diagonal 

matrix element equations again have a unique solution, but because T>£ is even there is now a condition on 
the solvability, with a periodic solution, of the diagonal ones: 

(2.36) iy (Vl 7q ,'y g )ds = f (q 1 + ^,...,q n + ^). 
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These conditions (2.36) determine f Q and hence the operator f (h, ■■■,I n )- To analyse the properties of f Q 
and of Qi, we use (1.1b), (2.16b) and (2.31c) to conclude that 

X>J = V + adQi{adQiD s + Vi) 

G C°°(R,£°)D 2 S + C°°(R,£ 2 )D S + C°°(R,£ e 4 ). 
We also observe that (2.32) is equivalent to 

(2.37) f (I u ...,/„) = J V t *vl\ V t dtds 

where T" = R"/Z™ is the n-torus, where t e T™ and where 

V(*i, ...,i„) := expitih expit n I n . 

Since V t belongs to the metaplectic representation, we have by metaplectic covariance of the Weyl calculus 
that 

(2.38) . f (h,...,I n )G V 2 c£ 4 
We then can solve for Qi m the form 

(2.39) Qi(s,y,D y ) :=-iL[J{T>$ - f (I u I n )}\ ds] e V I /0 (M 1 ) <g> £f 

where J is the indefinite integration in s with zero mean value and where the operator on the right is 
extended to L 2 (M} x R n ) using (2.35b). 

To make the integration more concrete and to analyze the solution, we again express everything in 
terms of complete Weyl symbols. Thus we rewrite (2.35b) in the form 

(2.40a) L~ 1 d s Qi(s, z, z) = -i{v\ | (s, z, z) - /o(ki| 2 , • • • , \z n \ 2 )} 

or equivalently 

(2.406) Q 1 {s,z,z)=Q 1 (0,z,z)-iL f [vl\ (u, z, z) - / (|zi| 2 , . . . , \z n \ 2 )]du 

Jo 

and solve simeltaneously for Qi and f a . The consistency condition determining a unique solution is that 

f L ± 

(2.41a) Q 1 (L,z,z) = Q 1 (0,z,z)-iL / [VS \ (u, z, z) - /o(|^x| 2 , . . . , \z n \ 2 )]du. 

Jo 

or in view of the twisted periodicity condition 

r L i 

(2.416) Q 1 (0,e ta z,e- ta z)-Q 1 (0,z,z) = -iL{ / VI \ Q (u, z, z)du - Lf (\ Zl \ 2 , . . . , \z n \ 2 )}. 

Jo 

As before we use that T>£ \ (u, z, z) is a polynomial of degree 4 to solve the equation. We put 

(2.42) Qi J2 1l;mn(s)z m z n , f (\ Zl \ 2 ,...,\ Zn \ 2 )= C ok\z 



. |2fc 
'ok\' c 

|m| + |n|<4 |fc|<2 



and 

I _ i _ _i i r L i 

T>S\o(s,Z,z)du := 2^ dg- mn {s)z m z n , dlmn-.^—j dS;mn(s)ds 

|m| + |n|<4 Jo 

As above, we can solve for the off-diagonal coefficients, 

(2.43a) g 1;m „(0) = -iL 2 (l - e '(— ")«)"! d% mn 

and must set the diagonal coefficients equal to zero. The coefficients c Q k are then determined by 

(2.436) c ok = d\. kk . 

It is evident that Qi and f (h, ■ ■ ■ , In) are even polynomial pscudodifferential operators of degree 4 in 
the variables (x, D x ), that Qi is weightless under metric rescalings and that the coefficients c k are of weight 
-2. They are essentially the same as the residual QBNF invariants. 
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We now proceed inductively to define self-adjoint polynomial pseudodifferential operators Q ±, fj(h, I n 

unitary /i-pseudodifferential operators W h i :— exp(ih^Qi) and approximate semi-classical normal forms 

such that: 

(2.44) 

i) T>h ■■= W*.W*^ . . . W* h V h W hh . . . W hh ^W hi ; 

ii) v\ ~ h- 2 + h~ l D s + E~ =3 h- 2 +%vl^; 



v: 



= [£>a,Q*-l]+2?_ 2 - 



(iii) for k > n — 2, 2? 2 2 _ 

(iv) **k>n-2,1>l 9 \ { =J^-^ l=Jj j + i 

(v) v\_^ G C°°(R,£™- 4 )£>2 + C°°(K,£ £ m_2 )^ + C°°(M,£ e m ); 

(vi) Q* G C°°(R,£ e fc+2 ); 

(vii) the conjugates under fj,(r a ) of the above operators are periodic of period L in s; 

(viii) fj(h, ...,/„) G 7>i +2 ; (is) «;«/*(/,•) = -2. 

To check the induction, let us assume these properties hold for k < N. Then 



(2.45) 2V =^-+l^ 2 ^-+l = 



m=2 .. - ^(adQ^) m (^) + E™ =1 E^3^ m ^ +f - 2 ^H0^±i) m (^|)- 



We note that 



[r 1 D s ,Q £±1 ]+^ 



and that further conjugations will not alter this (or lower order) term(s), proving (hi). 
When N + 1 is odd, Q w+i must solve: 

2 

(2.46a) {[L" 1 J D s ,Q«±i]+I?_ 2 ^i}|o = 

2 2 

while if TV + 1 is even Q w+i and f w-i must solve 

2 ' 2 

(2.466) {[i _1 I> a ,OiC±i] +^ f «^i}|o = ...,/„)• 

2 2 2 

JV 

If iV + 1 is odd, (2.45v) implies that the diagonal matrix elements of 2? 2 JV _ 1 | Q vanish. Hence, as in the 

2 

case of Qi, the solution of (2.46a) is given by 



(2.47a) 



Qi 



\o ■= ~iL J V 2 N _ 1 \ ds 



where the constant of integration Q w+i | o (0, z, z) is defined so that the solution satisfies the periodicity 
condition analogous to (2.28). In the even case, as in the case of Qi we set 



(2.476) 

or equivalcntly 
(2.47c) 

and 
(2.47d) 



1 1 1 f L n. 

fa_±(qi + -, ...,q n + -) := - / (V 2 N - t 7 g , J g )ds, 



fa^i (Ji,. ..,/„) := y [ L I V t *V^\ V t dtds, 

2 L Jo JT" 2 



Q 



n+i\ := —L j {D 2 n-i - fz=±(h, ...,I n )}\ ds. 
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As above, Q n±i \„ is then extended to all of TL a thru (2.31b). The indefinite integrations can be precisely 
defined, as above, by expressing everything as a polynomial in z, z and solving the resulting algebraic 

N 

equations for the constant of integration. By (2.47 c-d) the solution has the parity of T> 2 N _ 1 , i.e. the parity 
of N — 1, and by (2.45) the parity property propagates to the case N + 1. The formula (2.47c) also shows 
that for even N — 1, f w-i (h , ....In) has the same order and weight properties as V 2 JV _ 1 . Since the latter 

2 2 

may be written in the form aD 2 s + bD s + c, with c e C°°(R, £^ +3 ), the former lies in Vj n f^ 3 , implying 
(viii). The formula (2.47d) then implies (vi). The periodicity property (vii) is maintained throughout. 
It follows immediately that 

(2.48) 2 - 2 -2 

= [D„Q f _i]+X> 2 i,. 

Conjugating back under /i(r Q ), we see that statement (i) of the Lemma then follows from (2.48) and from 
(2.44iii-iv) and statement (ii) then follows from (2.44 viii). The weight property (ix) is visible from (2.47c). 

□ 

As a transitional step to the quantum Birkhoff normal form, let us rephrase the above Lemma in 
terms of the actual Fermi normal coordinates and inverse Planck constants {r kq } in place of (Lh)" 1 . We 
thus introduce the space 0*(N 7 ,T,{r k q}) of semi-classical Hermite distributions associated to the non- 
homogenenous isotropic manifold T := R + 7, i.e. 

T := {(s, a, y, rj) e T*(N 7 ) : a = 1, (y, V ) = (0, 0)}. 

It is, by definition, the union of the spaces m / 2 of elements of order ^,ra € N, given in Fermi normal 
coordinates by asymptotic sums of the form 

oo 

u(s,y,r^) ~ (r fc ,)*e ir *« s 5Z(rfc,)-*/„(a )V ^i/) 

n—o 

with /„ e C|°(K, 5(M")). Here, <S(M") is the Schwarz space, and the subscript q denotes the space of 
functions of this form satisfying 

f(s + L,y) = e iK « s f(s,ty). 

Aside from the restriction to Schwartz functions and the half- integral orders, it is the image under the 
rescaling operator T( r )-i of the space O* (iV 7 , A, {rk q }) of semi-classical Lagrangean distributions associated 
to the non-homogenenous Lagrangean 

A := graph(ds) = {(s, a, y, rf) : a = 1, T) = 0.}. 

Similar spaces of oscillatory functions could be, and implicitly are, defined in the model S* 1 x R" but the 
scaling aspect comes out most naturally on iV 7 . (For background on semi-classical Lagrangean distributions, 
see [CV2]). 

Let us denote by O* {N 1 , A, {n^}) the subspace of elements of O* (S 1 x R", A) satisfying Cu(s, y, r^ 1 ) <~ 
0. Let us also denote by Wk q the transfer of W -i above to the normal bundle, i.e. 

kq 

W kq := n(a s )W -ifj,(a s )*. 

kq 

For the sake of simplicity we will be a little negligent here of the rescalings. Under this operator, the kernel 
of C goes over to the space 0*(A 7 , A, {r^}) of elements annihilated by W£ q £Wk q - It is clear from the above 
lemma that this space is stable under Wk q * A r -iWk q . We may interpret what was proved in the lemma as 

kg 

giving a semi-classical normal form for A or A r -i in the following sense: 

(2.49) Theorem The r^ q -pseudodifferential operator 

W kq : 0* (N 7 ,A,{r kq }) -> 6* (N^ A, {r kq }) 

has the properties: 
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WW, 1 - / ~ Wkq W* kq -7-0 
^T-\ AT r _, W fe9 ~ £ 2 + / (7 7l , / 7 „) + /l(/7l ';-' V) + . . . . 
WV^^M = e irk <- s W kg U q (s,y)$ kq (s,y) 

(see (2.1)). 

3. Normal form: Proof of Theorem B 
The intertwining operators W kq will now be assembled into the Fourier-Hermite -series integral operator 

(3.1) W 7 : L 2 (Sl x E n , dsdy) -» L 2 (Sl x R", dsdy) 

^ 7 J] />,9)e^ s t/ 9 (s,?/)= ^ />, g ) e ^ s ^ fcg C/ 9 (s,y). 

(fc,?)6N"+ 1 (fc,g)£N"+ 1 

Also, the dilation operators will be assembled into the operator 

(3.2) T : L 2 (Sl x R™, dsdy) -» L 2 (5i x R™, dsdy) 

T ^ f(k,q)e ir - a U q ( S ,y) = ]T />, ^^(s, ^j/). 

(fc,g)GN™+ 1 (fc,g)£N»+ 1 

By theorem (2.49) we then have, at least formally, 

(3.3) W~ 1 T- 1 ^TW 1 ~ C 2 + / (7 7l) 7 7n ) + llihlipM +.... 

The purpose of this section is to make this equivalence precise. Since it is independent of the model we will 
carry out the proof in the basic model and use the weightless Fermi normal coordinates (s, x). For the sake of 
notational simplicity the W's and Q's, transferred back to the model in this way, will be written as Ws and 
Q's. Also in place of the 's we will use the weightless hu q s with = (2irk + Y^j=i (ij + \) a j) = Lr kq . 

(3.4) Propostion TW^T^ 1 is a (standard) Fourier integral operator, well-defined and invertible on the 
microlocal neighborhood (0.1) in T*{S\ x R"). 

Sketch of Proof: 

To simplify, we first consider the unitarily equivalent operator TWT^ 1 in the microlocal neighborhood 
(0.1) in the twisted model, with 

(3.5) W : l~t a — > H a 

W{e ll± ^ lq ) :=e lT ^W hkq j q , 

and with f the dilation operator like (3.2) relative to the basis e Lhk "^ q . We then factor TWT^ 1 as the 
product TWf- 1 =j*fvf~ 1 where: 

(3.6) V : l~t a — > Lf oc (R x R x R") 

V = nf =0 exp[iD- i Q i {s l ,y,D y )} 

that is, 

Ve lTi ^ lq {x) := e lTi ^W hhq {s\x,D x ) lq {x), 

and where 

(3.7) j* : C°°(R x R x R") -> C°°(R x R") 

ff{s,x) = f{s,s,x) 

is the pullback under the partial diagonal embedding. We note that the variable s' in V occurs essentially 
as a parameter, so we can (and often will) regard V as a one parameter family of operators V s i on TL a . 
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It is not clear that the infinite product in (3.6) is well-defined, nor what kind of operators the factors 
are. On the second point, we note that D s is elliptic in the set (0.1) in T*(S} j x 1") and that, as an 
operator- valued function of s', 

(3.8) TD-^Q^s'.x.D^f- 1 e C°° (Sl,* 1 ^ x R n )). 

Here of course C°° (S^ , -4) denotes the smooth functions of s' with values in the algebra A; to simplify 
the notation we do not indicate the possible twisting. Indeed, we first observe that conjugation by T 
converts the mixed polyhomogeneous- isotropic algebra ^*(S^)^)W* into the pure polyhomogeneous algebra 
V*{Sl) <g> **(R"). This essentially follows from the fact that 

Ta™(x, A^T- 1 = a w (D2x,D7*D x ) 
(see [G.2] or [CV] for further details). Obviously, TD.f' 1 e ^(Sl x MP), and by (2.44(vi)) we also have 

TQ 3 _f* g C°°(6'i,^' +1 (S' 1 x IP)). 

_ i ~ 

A simple calculation shows further that TD S 2 Q i(s' , x, D X )T~ is a one-parameter family of pseudodiffcr- 
ential operators of real principal type over 5^ x M™ whose principal symbols have nowhere radial Hamilton 
vector field. It follows in a well-known way that these operators are microlocally conjugate to D\ (see [Ho 
IV, Proposition 26.1.2]) and hence that their exponentials give a smooth one-parameter family of Fourier 
Integral operators. Consequently, after conjugation with T, the factors in (3.6) are for each s' G S 1 microlo- 
cally well-defined Fourier Integral operators in the neighborhood (0.1), with smooth dependence on s' . As for 

the infinite product, we note that the principal symbol of TD S 2 Qj_T~ 1 has the form cr~iqi(s', \foy, 
where q± is a homogenous polynomial in (y,rj) of degree j+2 (2.36(iv)). Its exponential can therefore be 
constructed microlocally in (0.1) as a one parameter family of Fourier integral operators over 5^ x 1" with 
phase functions of the form 

1 (s', s, s"a, x, x", £) = (s- s")a + (x- x", £) + Vi (s', s, s", a, x, x", £) 

2 2 

with 

(3.96) Vi e 0, +2 5 c 1 i 

Here, denotes the space of classical symbols of first order in (0.1) and Ok denotes the elements which 
vanish to order k along (x, £) = (0, 0). (Henceforth we will use the symbol Ok for objects of any kind which 
vanish to order k on this set.) It follows that for fixed s' and sufficiently small e, the phase parametrizes the 
graph of a homogeneous canonical transformation of T*{S\ x M") which is C j+1 - close to the identity. Any 
finite product in (3.7) is therefore a clean composition and the phase V&at of the product parametrizes 

the corresponding compostion \n of canonical transformations. By (3.9b), we have xn+i — Xn € Ow+i 
and so there exists a smooth one parameter family homogeneous canonical transformations \s' , equal to the 
identity along (a;, £) = (0,0) such that \s' — Xs'.n € On+i and with a generating function i/'oo satisfying 

ipca — 4>n € On (cf. [Sj, Proposition 1.3]). Regarding the amplitudes, the situation is of a similar kind. 

— i ~ 

Denoting the amplitude of Texp(D s 2 Qi)T~ by ai(s' , s,a,yr]) we have 
(3.10) ai = 1 mod0 3+2 S°. 

It follows that the amplitudes a,N of the products H-*L satisfy a^+i — a at € On+iS° and hence there exists 
an element a^, satisfying — G On+iS°. The pair (aoo, V'oo) then determine a one-parameter family of 
Fourier integral operators of order zero which agrees with n^L Q up to an error which also vanishes to order 
N. Finally, composing with j* just sets s' = s in the kernel, which visibly remains Fourier integral, with 
phase function parametrizing the graph of a canonical transformation \- Invertibility then follows from the 
fact that W*W is a pseudodifferential operator with principal symbol the square of j*(oy), which is easily 
seen to equal 1. 
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The proposition follows by expressing 

TW^T- 1 = Tfi( a yf- 1 fwf- 1 ffi(a)T- 1 

and noting that T^(a)T _1 is also a standard Fourier Integral operator. □ 
We now give the proof of the quantum normal form Theorem B for y/A, stated in an equivalent form in 
terms of W 7 . As in the introduction, the notation A = B means that the complete (Weyl) symbol of A — B 
vanishes to infinite order at 7. We also use the notation Oj^ m for pseudodifferential operators of order m 
whose Weyl symbols vanish to order j at (y,rj) = (0,0). Here, pseudodifferential operator could mean of 
the standard polyhomogeneous kind, or of the mixed polyhomogeneous-isotropic kind as in ^ k (Sj^) ® W l , 
in which case the total order is defined to be m = k + 1. To simplify notation, we will denote the space of 
mixed operators of order m by ty™ x {S\ x ^")- 

(3.11) Theorem B Let TW^T^ 1 be the Fourier Integral operator of Proposition (3.4), defined over a 
conic neighborhood of R + j in T*(S£ x R"). Then: 

W- 1 T- 1 yfKTW 1 = Pi(r,/ 7 i, J 7 „) + P (£,/ 7 i, I jn ) + -.., 

where 

(3.146) Pl {C, J 7l) / 7 „) ee C + P' 2 '(W-,V) + jfftW.V) + 



P— m (£ 1 I^fl j ■ • ■ 7 lyn ) — ^ ^ 



00 [k— m] 



(Ley 



k—m 

with p\ k m \ for m= -1,0,1,..., homogenous of degree l-m in the variables (7 7 i, 7 7 „) and of weight -1. 
Proof: 

As a semi-classical expansion in the "parameter" h = (3.3) may be rewritten in terms of \/A : 

(3.12) w-^Vatw^ ~ c + Pl(/71 ; 7 /7 " } + pa(/7i, ■-, A») + _ 

From the fact that the numerators /j(7 7 i, J 7 „) in (3.3) are polynomials of degree j+2 (2.36viii) and of 
weight -2, the numerators Pk(Ifi, I 7 „) are easily seen to be polynomials of degree k + 1 and of weight -1. 
Hence they may be expanded in homogeneous terms 

(3-13) P k =P l k k+1] +P l k k] + ... Pk °\ 

with pjj?' the term of degree j and still of weight -1. The right side of (3.12) can then be expressed as a sum 
of homogeneous operators: 

(3.14a) Px(£, J 7l , / 7 „) + P (C, J 7l , I 7 „) + . . . 
with 

/o i/n p(r t t ^— r 1 Pi ^7i' ^7") 1 P2 |(jrjj ^7») 

(3-146) -m.(1-, i 7 i, i 7 „) = L-\ — 1 — p h... 



00 [fc— m] / 7- t \ 

P f£ / , 7 ) = V ^ Khl,-,hn) 



k—m 

The remainder term in (3.12) can be described as follows: 
(3.15) W- 1 T- 1 yfATW 1 -[C + Elillli^iM + Mlju^n) + . . . + Pm (/ 7 i, J 7n ) ] 

e ©fcL + o 1 2(m+1 - fe )*^ fe (^ x M"). 

The remainder terms in (3.14b) are given by: 
(3.1ba) J 7 i, J 7 „) - [LH — 1 1 (LQk J 
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e 02(^+2)^(51 x E") 



(3.166) P_ ro (r,/ 7 i, V) - £ ^ y^;" 7 " j G 2(JV+ i- ro )*^(5i x E"). 

Hence the expansion (3.12) is also asymptotic in the sense of = . For the statement of Theorem B in 
the introduction, we only need to conjugate under /u(a). □ 
Remark The remainder in (3.15) could be equivalently described by saying that its complete symbol lies 
in the symbol class n S 1 ^" 1 ^ of Boutet-de-Monvel [BM]. 

4. Residues and Wave Invariants: Proof of Theorem C 



The aim now is to use the normal form of v A near 7 to calculate the wave invariants a 1 k associated 
to 7. In terms of the model (see the statement of Theorem B in the introduction) we may write the normal 
form as 

(4.1) V=^{h,...,I n )[V N +B N ] 
with 

~ „ , Pl(h,-,I n ) , , p N (h,-,In) 

VN - =1Z + UJZ + '" + (LRT 

and with 

Bn e ®^0 2(N+1 _ k) ^- x \S\ x R»). 

Our first observation is that we can drop a sufficiently high remainder term Bn in the calculation of a given 
wave invariant. We prove this by working out, roughly, which parts of the complete symbol of a general first 
order pseudodiffcrential operator P of real principal type contribute to the wave invariant of a given order 
associated to a non-degenerate closed bicharacteristic 7. 

In the following proposition, (s, y) will denote any coordinates in a tubular neighborhood U of the 
projection of 7 to M with the property that the equation for R + 7 in the associated symplectic coordinates 
(s,a,y,t]) is y = 77 = 0. Also, Op we will denote a fixed quantization of symbols in a conic neighborhood 
(0.1) of M+7 to pseudodifferential operators, and p(s, a, y, 77) ~ p\ +p + . . . will denote the complete symbol 
of P. The Taylor expansion of pj at M + 7 will be written in the form, 

Pj{s,<r,y,v) = v 3 Pj(s,l,y,-) = cr J (p' 01 +p\ 1] + ...) 

with p L J n \s 1 1, y, 2) homogeneous of degree m in (y, ^). We set Pj := Op(pj), Pj" 1 ' := Op{p^) and P^ N = 
T, m <N p j m] - Finally, we will denote by r 7fc (P) the coefficient of (t- L + iO) k log(t- L + iO) (or of (t-L + iO)- 1 
in the case of k — —1) in the singularity expansion of the microlocal unitary group e ttp near 7. 

(4.2) Proposition 

r 7fc (P) = r jk (P^ + P?*W + ... P° fe _!)- 

Proof: As mentioned several times, the wave invariant T 7 fc(P) is given by the the non-commutative residue 

(4.3) r 7fc (P) = resD k t e ltp \ t=L = resP k e ltp \ t = L 

(see [Z.l] or [G.2]). We first describe how (4.3) leads to a local formula for r 7 fe(P) in terms of the jets of the 
amplitudes and phases of a microlocal parametrix 

(4.4) F 7 (t,x,x')= [ e i ^ t ' x ' x ''^a{t,x,x',e)i){x,e)de 

for e ltp near {L} x R + j x R + j. The remainder of the proof will connect this data to the complete symbol 
of P. 

The amplitude a in (4.4) is a classical symbol, hence has the expansion a ~ Y^jLo a -i W1 tti a_j 
homogeneous of degree — j for |0| > 1. Since the residue or wave invariant depends only on the wave kernel 
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in a microlocal neighborhood of {L} x 7 x 7 it is unchanged by the insertion homogeneous cut-off ip, supported 
in (0.1) and equal to one on a smaller conic neighborhood of 7. The phase <p(t, x, x' ' , 9) parametrizes a piece 
ofthe graph of the Hamilton flow exptH pl of p\ near {L} x 7, and may be assumed to consist of only n phase 
variables. Writing x — (s, y) and 9 = (a, 77) as above and choosing ip of the form we get 

(4.5) r 7fe (P) = Res s=0 [ [ l\( k HL, S ,y,s,y,a,7 1 )^)e^ L ' s -y^y^\a\- s dadr 1 dv(s,y) 

Jr+ J-Rn-i Ju a 

where for certain universal constants C a p, 

(4.5.1) «(*>:= ]T C^(Df^)^ . . . (Df^-D^a. 

|a| + |/3|=fe 

We note here that the residue can be calculated using any gauging of the trace, in particular by powers |cr|~ s 
of the elliptic symbol a in (0.1). Changing variables to fj := % and denoting by o^-j the term of order k — j 
in the polyhomogenous expansion of the kth order amplitude a^ k \ and by B t the ball of radius e in M5, we 
get 

(4.6) r lk (P) = 
£ Res s=0 J J r+ J e ^.'.»^.fl) s W(L, s, y, s, y, 1, rj)i>(s, y, fj)<r n+k -^- s <dadfjdv(s, y) 



3 = 

OO 



E/ / Tpaf\<j>(L,s,y,s,y,l,ri) + iOy + i- n - k dfjdv(s,y). 



j=0 JUJB e 

The non-degeneracy of 7 as a closed geodesic implies its non-degeneracy as the critical manifold of <j>{s, y, s, y, 1, fj). 
Hence we have 

4>{L, s, y, s, y, 1, fj) = 4>(L, s, 0, s, 0, 1, 0) + Q s (y, fj) + g(s, y, fj) 

with Q a non-dengenerate symmetric bilinear form and with g = s (|(y, fj)\ 3 ). Expanding 

00 

(Q s (y, fi)+i0 + gf = C { x, P )(Qs(y, fj) + i0) x ~ p g p 

p—o 

we get that 

(4.7) r 7fe (P) = C {s+j _ n _ k , p) Res s = / / / (Q s (y,fj) + i0y +j - n - k - p (g p a k %i>J(s,y))dsdydfj 

j,p=o J~/ J\y\<e J B e 

with J(s, y) the volume density and with binomial coefficients C(A,p)- The family of distributions(Q s (?/, fj) + 
i0) x is meromorphic with simple poles at the points A = — (n — 1) — r, r = 0,1,2,... and with residue 
C r (Q~ 1 (D y , D rj )) r S(y, rj) for certain constants C r (see, [G.S., p. 276]). Here we have dropped the tilde in the 
notation for rj. For s = there are possible poles when k + p + 1 — j >0 with residues 

(4.8) C r C {j _ n _ kiP) f g7 1 (^,A ) ) fe+P+1 ^(3 P 4 fc ij)l(y=,=o)^. 

However, the residue vanishes unless 2(k + p + 1 — j) > 3p since g p vanishes to order 3p, constraining the 
sum in (4.6) to 2j + p < 2(k + 1). Hence we have, with new constants Cj kp , 

(4-9) r lk {P)= C ^ [ Q^\D y , D v ) k+p+1 -^g p a k %J)\ iyiV)H0fi) ds. 

j,p:2j+p<2(k+l) J ' y 

It follows first that j < k + 1, hence that only the terms a k k \ . . . , a^ k \ in the amplitude contribute 
to (4.9). To determine the corresponding range of terms cij in the amplitude a, we observe from (4.5.1) that 

(4.10) a k %:= Y, C a p{D^...{D^)^Dt^a k _ j _ m , 

\a\ + \f3\=k 

hence that only the terms a a , a_i, . . . , a_fe_i contribute to (4.9). 
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The same kind of calculation determines resAe ltP \t=L, for any fcth order pseudodifferential operator 
A supported microlocally in (0.1) in terms of the complete symbol of A and the phases and amplitudes in 
(4.3). In particular, we see that resD\ Ae ttp \ t= L = if A has order — (k + 2). This implies: 

(4.11) resD k e itp = resD k e it{Pl+ - p -^ ) 

since e itp = A k (t)e it( - Pl +-+ p -("+^ with A k (t) - I £ *- fc - 2 . Indeed, with H k = P 1 + . . .P_ {k+1) ,V k = 
P - H k ,B{t) := e- ltHk B{t)e ltHk , we have 

D t A k = V k {t)A k {t), A k (0) = I. 

Thus A k (t) is the time-ordered exponential of a pseudodifferential operator of order —k — 2, hence equal to 
/ modulo v]/- fe - 2 ; consequently so is A k (t). 

We next consider which jets j™a_i of a Q , . . . a_fc_i and j™<fi of 4> contribute to (4.9). Here by j™a_j 
we mean the Taylor polynomial of degree m of a_i(i, s, y, s, y, <r, -q) in the (y, 77) variables for t near L. From 
(4.9-10) it is evident that the maximum number of {y, rf) derivatives on a_i occurs when p = in (4.9) and 
when \/3\ = k,k — j — |/3| = —i in (4.10). It follows that at most 2(k + 1 — i) derivatives fall on a_i in 
(4.9). Also, the maximum number of (y,rj) derivatives on the phase occurs in terms where a factor of g is 
differentiated the maximum number of times; namely in terms with p = 1 and with j = 0, in which the 
phase is differentiated 2(k + 2) times. Hence 

(4.12) Tjk (P) depends only on ff k+ Va ,j* k a- 1 ,...j°a_ (k+1 y, % {k+2) <P- 

To give bounds on the jets of the terms pj(j — 1, 0, • • • — (k + 1)) in the symbol of P which contribute 
to the jets j^ k+1 *'a_i, we now have to consider some details of the construction of the parametrix (4.3). 
We first recall that the amplitudes are obtained by solving transport equations of the form 

j j-m 

(4.13) ^^[Pi 

m— u—o 

where V- Vt m = V-^m((j), s,7/, D s , D y ) is a differential operator of order < m obtained from Op(p- u ) as 
follows: 

00 

(4.14.1) e-MOpip-Je** ~ ^ p- v ~ m V- v , m 

<j> being the phase in the microfocal parametrix (see [Tr, Ch.VI,(5.28)]). Consistently with (4.11), only the 
terms pi,p , ■ ■ ■ ,P-( k +i) m the symbol of P contribute to the transport equation for a_j(i = o, . . . , — (k+ 1). 
To determine which jet ofp-^ is involved in "P— 1/. 771, we further recall [loc.cit.,Theorem 3.1] that the expansion 
in (4.14.1) is obtained by re-arranging terms in the expansion 

(4.14.2) e-^Op(p^)e^u ~ V -^ P - V {s, y, pd<j>)M a (cf>; p, D s ,D y )u 

a 

with 

N a {& p, D x )u = D«e^ (x - x 'K(x')\ x = x , 

(fo{x,rf) = <j>(x) - (j>{x') - (x - x',V(f>(x)}. 

In view of the fact that j\f a is a polynomial in p of degree we see that \a\ < 2m in any term of (14.4.2) 
contributing to V~ v ^ m - Hence V- V . m involves at most 2m derivatives of p~ v , and so the transport equation 
for a_i involves at most 2(i — v) derivatives of p_ v (y = — 1, 0, . . . , i.) 

To draw these kinds of conclusions about the amplitudes a-i themselves, as opposed to the coefficients 
in their transport equations, we have to take into account the initial conditions in the transport equations. 
Unfortunately, the defining initial conditions occur at t = while we are interested in the long time behaviour 
at t = L. Were there no conjugate pairs along geodesies near 7, the transport equations could be solved up 
to t = L and since the initial conditions can be taken in the form a \t=o = 1, &-i|t=o = we could conclude 
that the a_i's are integrals of polynomials in at most 2{i — v) derivatives of p- v and in at most 2(i + 1) 
derivatives of the phase (coming from the A/^'s). However, in the case of elliptic closed geodesies, there will 
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always be conjugate pairs for L sufficiently large and we cannot construct the parametrix so simply. Rather 
we will use the group property U(L) = U(L/N) N with N sufficiently large that a parametrix for U(L/N) 
can be constructed by the geometric optics method. We then have to eliminate all but n phase variables in 
the power, which will lead to crude but serviceable bounds on the order of the jets. 

We therefore begin with the construction of a short time parametrix Fgk, valid for \t\ < S, of the form 
(4.4) with phase <p(t,x,x',0) = S(t,x,9) — {x',6), satisfying d t S + pi(x, d x S) = 0,S\ t=o = (x,6), and with 
amplitude a <~ Y^j=o a -o satisfying (4.13) on \t\ < S and with initial conditions a o | t=0 = 1, a-j\ t=0 = 0(j > 
0). By the observations above, we have that ag t -i(t, x, x' , 9) involves at most 2{i — v) derivatives of p~ v {y < i) 
for t in this interval. As is well-known, we then have e ttp — Fsk(t) = R k(t) with the remainder Rsk(t) a 
Fourier integral operator of order — (k + 2) associated to the graph of exptH pl for \t\ < S. Since Fst is of 
order zero, it follows that for any TV and for \t\ < 5, e ltNP = Fgf.(t) modulo Fourier integral operators of 
order —k — 2 in this class. Hence, the desired parametrix (4.4) can be constructed by choosing N so that 
L/N < S and by eliminating all but n phase variables in Fgf. by the stationary phase method. Just as in 
(4.11) the remainder terms with a factor of R$k will not contribute to r 7 /j(P), nor will terms a_j in the final 
amplitude with i > k + 1. 

To complete the proof of the proposition, we have to count the number of (y, ^-derivatives of p- u which 
appear in the terms 

(4-15.1) A fc " l/ "Q7 1 (^.^) fc+1 " i afc-i-|/3|l(w,,)=(o,o). 

and the number of derivatives of the final phase (hence of p\) that occur in the terms 

(4.15.2) Q; 1 ^, AO^-Wr^ 1 • • • (A a ^) /3 "A a " +1 ^-| /3 ||(,,,)= ( o,o) 

with d-i the terms in the final amplitude. We start from the fact that, in an obvious notation, 

(4.16) F fe ^= / / e H<h+-M V a- il ...a- iN d0 1 ...d0 N dv(x 1 )...dv{x N - 1 ) 

and that the elimination of all but n phase variables by the stationary phase method gives essentially the 
same formulae for the final amplitudes d-i as in the case of non-homogeneous Lagrangean distributions of 
the form J u{x,y)e llAl ^ x ^dx with f(x,y) a real-valued function defined near (0,0), with f' x (0, 0) = 0, and 
with fx X (0,Q) non-singular. In this situation, we have (see [Hoi, Theorems 7.7.5-6]) 



/ 



u(x, y y uf< - x > y) dx = C V Lf j v u(x(y), y)uJ- ] 

\det(u,f» x (x(y), y ))\-i j+ f <M L °' V 



modulo terms of order 0(w ^ M ^), with Lfj tV a differential operator of order 2j whose coefficients involve at 
2j derivatives of /". It follows that d-i has the form (with I = (ii, . . . ,ijv), $ = <t>i + ■ ■ ■ 4>n)) 

^ ' LQ jCi—i 1 . . . &—i N . 

I,j:\I\+j=i 

Since also a-i has the form 

FiipuDpu . . . , D 2 ^ Pl , Po , . . . , D 2i Po , . . . 0, . . . , D 2i 4>) 

with Fi a multiple integral of polynomials, and with D k some differential operators of degree k, we see that 
also d-i has the form 

Kipu D 2 ^ +1 \p , . . . , D 2i Po , D 2 p- t+ll p- t - 0, . . . , £ 2 (* +1 V). 

Hence as regards number of derivatives of p\, . . . ,p~i, (j), the a_i's behave exactly as the a,'s, that is, involve 
at most 2(i — v) derivatives of p- v . Hence (4.9)-(4.10) apply to the final amplitudes, and we conclude that 

2(fe+2)„ ,-2(fc+l) -o -2(fe+2) . 



r 7fc (P) depends only on j^ k+2) Pi, j 2( * +1> Po, ■ ■ ■ j°P-(k+iy, j. 



Since j 2 ( fc + 2 )^ depends only on j 2 ^ k+1 ^pi, the proof of the proposition is complete.. □ 
(4.16) Corollary 

T 7 fe(X>) = r 7 fe(Dfc+i). 
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The following Lemma combines the previous results in a form applicable to the calculation of the 
wave invariants. The notation TrAP~ s (with P £ $' elliptic) is short for the zeta function obtained by 
meromorphic continuation of the trace from Res » to C. 

(4.17) Lemma 

r 7fc (VA) = Res s = D*Tr${n,h, . . . , I n )e iW "+ l TZ- s \ t=L 



Proof Since 1Z is elliptic in the essential support of ip, the trace on the right side is well defined and 
has a meromorphic continuation to C (cf. [G.2], [Z.1,5]). The residue is of course T jk (T>) by the previous 
proposition. Hence it suffices to show that 

(4.18) r 7fe (\/A) =r lk (V). 

This however follows from the previous proposition combined with Theorem B: Indeed, the proposition shows 
that T lk {A) = T 7/ t(_B) if A = B in the sense of Theorem B. Also, r lk {W AW\) — T lk (A) if W\ is a parametrix 
for W on the essential support of A, as follows from the tracial property resW AW\ — resW{W A of the 
residue (see [Z.1,5] for instance). □ 

We come now to the calculation of the wave invariants as residue traces of the normal form wave group. 
But we will simplify (4.17) further before evaluating the residue trace. First, we rewrite everything in terms 
of D s , and H a using (0.3-4). Since 

p v {Ii,...,I n ) p v {I\, . . . ,I n ) H a 1 H a 2 

(L1Z) V = (LD s y V- V LD S + 2 V{V - 1){ LD^ 

and since we can drop the D s ( fe+1 ) +I/ H^ +1 ~ u and higher terms by Proposition (4.2), V k+ i can be written 
in the form 

(A ^Q\ t> - n mtt i Pi(h,---,I n ,L) p2(Ii,...,I n ,L) p k+1 (h,...,I n ,L) 

(4.19) V k+1 =D S +H a + — + ^2 + ••• + 

modulo terms which make no contribution to t^ 7 . Secondly, we can use LD S rather than 1Z as the gauging 
elliptic operator in (4.17). To simplify the notation we will denote all but the first two terms on the right 
side of (4.19) by 

(A 0(X\ v tn t t t\. Pi(h,---,I n ,L) p 2 (h, ■ ■ ■ , I n , L) p k+1 (h,..., I n ,L) 

(4.20) V k+1 (D S , h, ...,/„, L) .= — + ^2 + • • • + IZD^+l • 

Then we have: 

(4.21) r 7fc (VA) = Res z=0 TrD^(D s ,I u . . . , I n )e u ^ LD ° +H ^ +v ^\LD s )-*\ t=L . 
We can now give: 

Proof of Theorem C: By (4.21) and the fact that e 2mLDs = I on L 2 (Sl) we get 

afe 7 = T 7 fc(\/A) = 

(4.22) Res z=0 Tr[y2(wLD s + H a ) + V k+1 } k e^e lLVk +^(^-)(LD s )- z 
with 

From the well-known spectra of D s , I\, . . . , I n we can rewrite (4.22) as 



(4.23) Res z=0 f2m z { ^ [i(27rm + ^{q, + ^)a 3 ) + V k+1 (m, Ql + 1 . . . , q n + ^)] h 

m=l (9i,...,gn)eN™ j=l 

. e iiP fc+1 (m, 9l + i,..., 9n + i,L)^_^) e i Ej = ife + iK}_ 



em 
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Regarding 

(4.24) V ^JLyTU^+h)^ 

as a smooth function ofagl",we can further rewrite (4.23) as 

n 

Res z=0 m- z {[-(2irm + ^ a 3 D » 3 ) + Pk+i (m, D ai ,...,D an , L)f- 

rn — 1 j — 1 



(4.25) e iLT k+1 (m,D ai ,...,D an ,L) £ ~r { J_ y ^ j=1 te + 5 )^ } 

(«i,.,fc)£N« £m 

Since ^{-^-) is for each to a finitely supported function of q, we can also rewrite (4.24) as 



x 



hm n" ^(^) Y e'E^i^to+i)"* 
a^o J - i v em ; ^ 

(4.26) = limn" ^(^)n™ _ _ . 

Recalling the definition of T(a) (0.12a), combining (4.23)-(4.26), and taking the limit as S — > 0, we get: 

1 n 

a lk = Res z=0 ^ m~ z {[—(2iTm + y~]ajD ai ) + Pk+i(m, D ai , . . . ,D an ,L)] k - 

rn—l j — 1 

(4.27) e iLV k+1 (m,D ai ,...,D an ,L)j 1 n T/_2tW( a )l 

J 6TO 

Here, we use the fact that T(a) is a a temperate distribution on R™ with singular support on U" =1 R x R • • 
2nZ • • • x R" (the factor of Z occurring in the jth position) and that the limit S —* in (4.26) can be taken 
in 5'. Since the cut-off smooths out the singularity, the right side of (4.27) is a smooth function of a and 
can be evaluated at the special values of a determined by 7. 

Since Vk+i(m, D ai , . . . , D an ,L) is a symbol of order —1 in to with coefficients given by polynomials in 
the operators D aj , we can expand the kth power in (4.27) as an operator-valued polyhomogeneous function 
of to. At least formally, we can also expand the exponential e tLV k+i{m.D ai ,...,D an ,L) j n a p 0wer ser i es an d 
then expand each term in the power series as a polynomial in mT 1 . Collecting powers of m, the right side 
of (4.27) can be put, at least formally, in the form 

(4.28) Re Sz=0 £ £ (£>«, L)?(^)T(a), 

m— 1 j—o 

with !Fk,k-j(Da, L) the coefficient of m i in (4.27). To justify this manipulation, we have to deal with the 
remainder term in the Taylor expansion of the exponential. We have 

e lx = ejy(ix) + rjy(ix) 

, , (ix) N 
e N (ix) = l + ix + --- + ±j^- 

rjy(ix) = {ix) N+1 b^{ix) 



b N {ix) = f ... f t%t^Z 1 1 ---t° e tNtN - 1 - t " ix dt -..dt N . 

Jo Jo 



Hence we have 



e iLV *^ := e N {iLV k+1 ) + (iLV k+1 ) N+1 b N {tLV k+1 ) 
with Vk+i short for Vk+i{m, D ai , . . . ,D an ). The ejy term of course contributes a finite number of terms 
of the desired form (4.28). For the remainder, we expand (iUPk+i) N+1 as a polynomial in mT 1 with 
coefficients given by operators QNp(D ai , . . . , D an ) and observe that each term has a factor of to _w_1 . For 
each such term, we remove the coefficient operator Qnp from the sum ^2 m , leaving only the factor of b^. 
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We then rewrite the resulting sum as a double sum ^2 mq as in (4.23), replacing all operators in D aj by their 
eigenvalues. Since b^iix) is a bounded function and since each term of the resulting sum has at least the 
factor m~ z ~ N ~ 1+ (possibly multiplied by a negative power of m), we see that the remainder is a sum of 
terms of the form 

(4.29) Res z = Q Np (D ai ,. . . , D a J V m -' +k - N - 1 - l b N (iP k+1 (m, q + )e^+^ a l 

kq 

We then observe that the sum is bounded by Ylm=i m~ Rez ~ N ~ 1+k+n , hence converges absolutely and uni- 
formly for Rez > — N + k + n. It follows that for N > (n + k) the sum in (4.29) defines a holomorphic function 
of z in a half-plane containing z — and since the operations of taking the residue in z and derivatives in a 
commute, each term (4.29) is zero. This justifies (4.28) and shows that it is actually a finite sum in j, say 
j < M (in fact M=(k+l)(n+k+l)). 

The residue in (4.28) is therefore well-defined and independent of e. Since ip(^)T(a) — > T(a) as 
e — > oo we must have 

oo M 

(4.30) a 7fc = Res z = Q ]T ]T m-*+ k ->F ktk _ j (D a , L)T(a) 

m—l j — o 

M 

= Res z = C(* + J - k)fk,k-j {D a , L)T{a). 

3=0 

Here, £ is the Riemann zeta-function, which has only a simple pole at s = 1 with reside equal to one. It 
follows that the only term contributing to (4.29) is that with j = k + 1 and hence we have 

(4.30), < hk =F k ,- 1 (D a ,L)T(a) 

completing the proof of Theorem C. □ 

5. Local formulae for the residues: Proof of Theorem A 

The characterization of the wave invariants in Theorem A is reminiscent of that of the heat invariants 
in [ABP] or [Gi], but involves non-local metric invariants near 7. We begin by determining the metric data 
which contribute to a lk . 

As always, we assume that 7 denotes a primitive closed geodesic, and denote the £th iterate of 7 by 
j e . As above, we use the exponential map along 7V 7 ~ Sj^ x R" to pull back the metric to a metric on a 
neighborhood of S\ x {0}, or more simply S 1 , in S\ x R n with the same wave invariants as g along 7. This 
reduces the theorem to the case M = S\ xl°, with S 1 a closed geodesic of the metric. 

We let J^l denote the manifold of m-jets along S* 1 of Ricmannian metrics on S 1 x M™ with the property 
that S 1 is a closed geodesic. We also let J™i ell denote the open submanifold of m-jets of metrics for which 
S 1 is non-degenerate elliptic. We also write J™ zi when we wish to identify 7 and 5 1 . The density I 7 k(s)ds 
of the fcth wave invariant is then a map 

where f2 1 (5 1 ) is the space of smooth densities along S 1 and where the jet order nik will be shown below to 
be mk = 2fc + 4. Since 7 7 fc(s, g)ds is independent of the choice of coordinates on S 1 x K" we may express it 
in terms of Fermi normal coordinates (s,y) with respect to a fixed normal frame e(s) for g. As usual, the 
metric coefficients gij(i,j = o, 1, . . . , n) will be understood relative to these coordinates. 

§5.1 The metric data in I kl 

We now claim that Jfc 7 is an invariant polynomial in the following data: 

(i) the curvature tensor R and its covariant derivatives V m i? with m < 2k + 2, contracted with respect 
to the Fermi normal vector fields and X^f=i c iaf~ with cj G R; 

(ii) the components Reyij(s), Imyij(s) of the normalized eigenvectors Yj e ® C of P 1 relative to 
e\, . . . , e n , and their first derivatives; 
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(iii) at most 2k + 1 indefinite integrals over S\ of (i)-(ii). 

(iv) the length L and inverse length L _1 of 7^; 

(v) the Floquet invariants (3j = (1 — e lQy ) _1 ; 

Indeed, by Theorem C, I 7 .k(s,g)ds is a density depending only on the data contained in T(a),L and 
in the coefficients of Tk.-i- The latter is identical to the data in the coefficients in pi, . . . ,pk+i, hence to 
that in p\, . . . ,pt+i- By Proposition (4.2), these depends only on jy k+4 g, and on the coefficients of the 
symplectic matrices r a (s), a s which arise in the intertwining of Ah to its normal form. The coefficients of r a 
depend only on L and the aj 's and those of a s depend only on the Jacobi field components , and their first 
derivatives yij . The second and higher derivatives of yij can of course be eliminated by means of the Jacobi 
equation. Regarding (ii), we recall (see [Gr, Theorem 9.21]) that the Taylor series expansions at y = of the 
9ij( s ,yY s , °f the co-metric coefficients g l; >(s,y), and of the volume densities and their powers, involve only 
the curvature tensor R and its covariant derivatives Vi?, . . . , V m R contracted with respect to the normal 
Fermi vector fields. Since the coefficients £2-^ of the semi-classical expansion of A/j depend only on this 
data, and since the intertwinings only introduce coefficients in data (ii)-(iv), the pfe(/)'s can only involve the 
gij's thru the curvature and its covariant derivatives. 

For the fact that I 7 fc(s) is a polynomial in (i)-(iv), it suffices to reconsider the construction of the normal 
form. Obviously the coefficients of the £'s and their metaplectic conjugates the 2?j's are polynomials in 
(i)-(iv). It then follows from (2.39c) that the fj(I)'s, hence the Pj(7)'s, are also polynomials in this data. 
Indeed, we argue inductively from the construction of the normal form that the coefficients change in the 

m m+1 

step from V£ to V L 2 as a result of taking commutators with operators whose coefficients are polynomials 

2 2 

T71+1 

in (i)-(iv). It follows that the coefficients in T> i : 2 are also polynomials in this data, with possibly two more 
-D s -derivatives due to the commutators with D s . Since the Jacobi data yij , i/ij enters in thru a linear change 
of variables, the degree of the polynomial in this data will be closely related to the degree in the (y,r]) 
variables, which is twice the degree in the I variables. The degree of the polynomial in the Jacobi data is 
however not necessarily the same as that (namely, k+2) in the /-variables: from the proof of Lemma (2.22), 
we see that the algorithm for computing the polynomials fj(Ii, ...,/„) involves taking operator commutators 
(or Poisson brackets of symbols); this lowers the order in the /-variables but not in the Jacobi data which 
are coefficients of the polynomials in the /'s. We will show below (see 'Jacobi degrees') that the order in the 
Jacobi data is no more than 6fe + 6. 

§5.2 Weights of 1^ and of the data 

To determine the weights of the polynomials in the data (i)-(iv), we now extend (and in part recall) the 
table in §1.4 describing how the various data transform under the metric rescaling g — > g e := e 2 g. As above, 
(s, y) always refer to Fermi normal coordinates relative to g, and for notational simplicity we put s = y a . The 
notations V, R, A refer respectively to the Riemannian connection, curvature and Laplacian. We distinguish 
A from the local expression Ym,j=o ®vi9 % * ^/9®Vj f° r ^ relative to the Fermi normal coordinate frame. 



9 


2 

eg 


9tj =9(9 Vi ,d yi ) 


9H = e 2 9(e~ L d yi ,e~ L d yi ) 


DZv9ij 




V, R, A 


V, R, e-' z A 


— >^;:, ..<>:, n'\9<>,, 


< J 7-7>T., ./>:< 9'\Jl<> :/ 


L, (s,y) 


eL, (es,ey) 


Vij ■= 9{Yi,ej) 




Kij = g(R(d s ,ei)d s ,ej) 


e-*Kij 



The trace of the wave group thus scales as 
(5.2.1) Tre lt ^ -> Tre 1 ^ 
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from which it follows that 

(5.2.2a) a 7fc = resVA^^ -> e~ k a jk . 

This can also be seen from the fact that a lk is the coefficient of (t — L + iO) k log(t — L + iO) which is 
homogeneous of degree k modulo smooth functions of t. Since a 7k = I lk (s)ds, and the integral scales like 
ds, we also have 

(5.2.26) I 7k (s,e 2 g) = e- k+1 I lk ( S ). 

As a check on the normal form, let us verify that wgt(a-y k ) — — k directly from a weight analysis of the 
normal form. It is obvious that 1Z = j^(LD s + H a ) scales like e _1 and so does each term in the expansion 

VA~n + Mh --> In) +.... 

Moreover the Pk{h, ■ ■ ■ , AO's and TVs are also of weight -1, as determined in Lemma (2.22) and Theorem 
B. Since the transition from p v to the p v +rS only involves multiplication of p v with the rth power of -^g- it 
is obvious that p v also has weight -1. Expanding the exponent in the residue calculation, we get that 

fe+l .jv -I 

(5.2.3) a lk = m™^ ' + H ^ + n} k (LV k ) N e^ 

with [j-(LD s + H a ) + V k \ of weight -1 and with (LV k ) of weight 0. Hence wgt(a 7k ) = —k. 

§5.3 Wave invariants and QBNF coefficients 

They are related as follows: The terms in (5.2.3) contributing nontrivially to a lk have the form 
(5.3.1a) L- fc +^+'--^i res ( i ^ s )^-(ji+2j2+...(/ s +i)j fc+1 )^ pjipj2 _ _ .pj*+i e «r- 

with 

(5.3.16) j a + j k+ i +£ = N + k, h + 2j 2 + ...(* + l)j k+1 =t+l. 

By the argument in the proof of Theorem C ( §4) , taking the residue removes the factors of LD S and replaces 
an expression res(LD)^ 1 F(I)e lHa by the value of F(D a + |)T(a) at a regular point. It follows that 

(5-3.2) a 7fc = J2 C*^£- fc+J ' 1+ - J '* +1 [HZ ftp* . . .ft£](Da + \)T{a) 

where the sum is taken over the indices specificed in (5.3.1b). 
§5.4 Jacobi degrees 

We now show that the degree of I lk as a polynomial in the Jacobi data yij,yij is < 6k + 6. The proof a 
detailed review of the construction of the polynomials fj{I\ 1 • • ■ , I n ) in Lemma (2.22) as well as the relation 
between the Jacobi degrees of the polynomials fj's and of the operators Tk-\ in Theorem C. For the sake 
of brevity, and since it is quite routine, we will be a little sketchy in a few of the details. We will use the 
notation J — deg for the degree of a polynomial in the metric data above with respect to the Jacobi field 
components. 

(5.4.1) Lemma J — degpj + i = J — deg pj +i = J — deg fj = 6j + 6. 

Proof From the formulae (4.19) relating the Pj{h, ■ ■ ■ , irt)' s an d the Pj(h, ■ ■ ■ , In), and from the fact that 
yf- has Jacobi degree 0, we see that 

(5.4.2a) J — degpj = max {J — degpi, . . . , J — degpj} 

On the other hand, the relation between the p/s and fj's is given by comparing (3.3) and (3.12): 

2 x , 

(5-4.26) fj = -p j+1 + 2^ Pi+iPr+i 

i+r+2=;+l 

Let us assume J — degfj = 6j + 6: an easy induction using (5.4.2b) then shows that J — degpj+i = 6j + 6, 
and another using (5.4.2a) shows that J — degpj + i = 6j + 6. 
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Hence we must prove (5.4.2b). By (2.39c) we have that 

(5.4.3) J - degfj = J - degV 3 ^ . 

3+- 

To determine the latter, we need to recall that T>_^ is the h 3 -term in the expression 

Ad(e ihj+hQ ^i)Ad(e ihi ^) . ..Ad^ Q h)V h 

with Ad{V)A := V~ 1 AV. Using that Ad(e lB ) — e lad ( B \ and expanding everything in a formal h-sehes, we 
find that 2? 2 is the ft,- 7 -term in the series 

(5.4.4) J2 h%~ 2 (W + ifi+i{V) N i . . .{h^f? ad{Q J+ i_) N >+? . . . ad(Q ifiV 2 _~ . 



n,Ni ,...,JV 2 ,_i 
7 



We claim 



(5.4.5) Claim: Suppose J — degQm = 3m for m < 2j + 1. Then: J — degT>_^ 2 = 6j + 6. 

Proof: The h? term in (5.4.4) is the sum of terms with indices satisfying ^ — 2+Nj + ± (j+|) + - ■ -+Nx(^) = j- 
Multiplying by 6 and using the hypothesis we find that 

i 2j+l 

J - degVi+i = n + 6 ^ y7V f = n + 6[j + 2 - |]. 

m— 1 

The claim now follows from the fact that n > 3 in any term in the sum. □ 

j_ 

(5.4.6) Claim: Suppose J — degV 2 } 1 = 3j + 3. Then: J — degQ j+i = 3j + 3. 

~2 + 2 2 

Proof: Follows from (2.39c,d) which implies that 

(5.4.7) J - degQ , i = J - degV 2 _ L 

"2" 2 2 + 2 

□ 

We then show: 

(5.4.8) Claim: J - degQ L = 3j. 

2 

Proof: We prove this by induction on j. For j = 1 it holds by explicit calculation: from (2.27a), degQi = 

i 

degVi — degCi = 3. It follows then by Claim (5.4.5) (with j=0) that J — degT>£ = 6, and then by Claim 
(5.4.6) (with j=l) that J — degQi = 6. The rest of the induction proceeds similarly and the details are left 
to the reader. □ 

The proof of Lemma (5.4.1) is completed by combining (5.4.3) and Claims (5. .4.5) and (5.4.6). □ 
Finally we have 

(5.4.9) Lemma: J — dega^ < 6fc + 6. 

Proof: As in §5.3, the kth residual terms in [(D s + H a ) + T , k ] k e lHa e lLVk involve only the products 
of the p/s with 

(5.4.10) 3o+h + ...jk+i+t = N + k, j 1 +2j 2 + ...(k + l)j k+1 =£+l, £<k,N<k + l. 

Using that J — degp m < 6m, the Jacobi degree of any such term cannot exceed 6(ji + • • • + jk+i)- We claim 
that the maximum possible value, subject to the constraints (5.4.10), is 6k+6. Indeed, since ji + • • • +jk+i = 
I + 1 — (j2 + 2.73 + . . . kjk+i) < k + 1 — (j2 + 2.?3 + . . . kjk+i) the maximum value is achieved when j\ = k + 1 
and all other j m 's are zero. □ 

This completes the proof of Theorem A. □ 
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6. Quantum Birkhoff normal form coefficients 



In this section we give a brief summary of the algorithm for calculating the quantum Birkhoff normal 
form coefficients B^.^j- From these coefficients one could determine the wave invariants a T fc as described 
in §5.3, but since the B^-kj's are simpler spectral invariants we choose to concentrate on them. In the next 
section, we apply the algorithm to the calculation of the coefficients £ 7;o j m dimension 2. 

Preliminaries 

(6.1) Definition The quantum Birkhoff normal form (QBNF) coefficients are the coefficients of the mono- 
mials P := If 1 . . . J;( n = \zi\ 2 ^ . . . \z n \ 2 i n in the complete Weyl symbol pk(\z\\, • • • , \z\n) of the coefficient 
operator p~k(h, ■ ■ ■ ,I n ) of the model normal form of Theorem B: 

p k (\z\\, . . . ,\z\ 2 n ) = B r,k,M 2j - 

j£N":|j|<fc+l 

There is of course something arbitrary in the emphasis on the pfc(ii, . . . , 7„)'s here. The coefficients 
of the monomials of the complete symbol of the pk{P, ■ ■ ■ , I n )'s are equally spectral invariants and we only 
prefer the Pk(h, ■ ■ ■ , ^n)'s to maintain consistency with the terminology of [G.l-2]. Moreover, the coefficients 
of the operator monomials P in either the p^s or p's are also spectral invariants and in view of the relation 
between wave invariants and the QBNF (§4) , it is the operator coefficients which are most simply related to 
the wave invariants. 

The crucial point for the effective calculation of the QBNF coefficients is that they can be obtained 
from the coefficients of the complete symbol of the semi-classical normal form (SCNF) 

oo 

W* h TZ h W h \ -h- 2 +J2 h k f k {h,- • •,/„) 

3=0 

of Lemma (2.22) (after restriction to weight 0.) Indeed, after substituting h = 7£ _1 , the square of the QBNF 
is the SCNF. Hence, the key invariants are really the coefficients of the monomials in the complete symbols 

/fc(M 2 , . . . , \z n \ 2 ) = c i;kj\ z \ 2j 
jen^,\j\<k+2 

of the coefficient operators fk(h, ■ ■ ■ , I n )- In the remainder of the section, we drop the subscript 7 from the 
notation. 

On the operator level, the relation between the pkS, pus and f' k s is very simple: the operators 
Pk(h, ■ ■ ■ , In) of Theorem B are related to the operators fk(h, ■ ■ ■ , In) of Lemma (2.22) by: 

(6.2a) f k = jPk+i + Pip i- 

i+j—k.i,j>l 

while the operators pj are related to the operators pk by: 

Vl = C k,(h,...,j k ) Hl J l Pk 

k>l,Uu-Jk)-k+\j\=e 

for certain universal (multinomial) constants O n the symbol level, the relation is a little more 

complicated since one has to compose the symbols in the Weyl calculus. In the following we summarise the 
steps involved in calculating the coefficients of the monomials \z\ 2 ^ in the complete symbols of the /fe's. To 
distinguish on operator from its complete symbol we use the notation f(I\, . . . ,I n ) for the Weyl operator 
with complete symbol f(\zi\ 2 ,...,\z n \ 2 ). We also recall that the notation | Q refers to the weight part of an 
operator relative to D s : (A 2 (s, x, D X )D 2 + A\(s, x, D X )D S + A Q (s, x, D x ))\ a = A Q (s,y, D y ). We will assume 
that the C's and 2?'s have been expressed in terms of the weightless Fermi normal coordinats (s,x) of §1.5. 
The complex coordinates z are given by z = x + i£. 

Summary of the algorithm 
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Step 


Relevant data 


Step 1 


Express /w-i in terms of T>2-^'s 

9. 2 


Formula for f w-i 


L 

/n-i (7i, ...,/„) := \ J o J T „ V t *T> 2 N _ 1 \ V t dtds (diagonal part; N odd) 


Equivalently 


- 1 * 

/«-l(<J1 + 3, ••.,?„ + 2) : = T./o (V * N-l Jq,7q)ds 
2 2 


Weyl symbol 


/«_l(|z 1 | 2 ,...,|z„| 2 ) = iJ2 lkl< N+3CN^ .\Z\*3 


JV + l 

Recursion for I? 2 N _ 1 

2 


~ 2?E N-l ~ 2 ^ 

T> 2 N _ 1 = term of order h 2 in 

^ AT+2 h -2+^ V f_^ h ^l {[L -l Dsj Q^ l]+ } 

+ h-*+Sv}_ 9 + YZ=2 h m ^Q(adQ^y-(L^D s ) 
+ E^ 1 E™ 3 h m ^- 2 £iadQ^r(vf_r) 


Operator Transport : 
equations for Q's 


N 

iV+lodd: {[L- 1 D s ,Qn +1 ]+V 2 n _ 1 }\ o = 

N + l even: {[L" 1 ^, Q«+i] + Z> V^}| = I n ). 


Symbol 
Transport equations 


~ ? ' N_ 

N+l odd: L~ i d s Q n+i (s, z, z) = —iD 2 JV _ 1 \ a (s, zz) 

iV + 1 even: Z, 'ftOn+i (a, 0, 2;) = -i(XJ> 2 N _i }| G (s, 2, 0) + ifw-i (\zi | 2 , . . . , \zi\ 2 ). 

2 9 2 


Homological 
equations 


N + 1 odd: Qn+i (0, e ia z, e~ la z) - Q«+i (0, z, z) = -iL V^ N _ X \ Q (s, zz)ds 
N+l even: Q n+i (0, e la z, e~ la z) - Qn±i(0,z,z) = -iL (V^ N _ 1 }\ {s,z,z)ds 


Solutions: 


in 

~ ~ s — 

N+l odd: Qn+i (s, z, z) — Q n+i (0, z, z) — iL J o V 2 N _ 1 (£, z, z)ds 

~ 2 ~ 2 2 

N+l even: Q n+i (s, z, 2) = Q n+i (0, z, z) 

w 2 2 
- iL fo{ V 2 N-i\o{t,z,z) - /iv-i (|zi| 2 , |z„| 2 )}| ds 


Step 2 


Express T>2-^ 's in terms of metric data 


Conjugate to £ 2 -f 's 


P 2 -t =M(^£*)£2-tM^*)" i 


r 

£ 2-f - 

_2 -term of 


( V, T \—'2, ~oo 1 o^V^ r\-l n oo a 1 , ■ /' r \ -lro 

+^ _1 (Er ? -i 9}Ld Xi d Xj ) + h-i(Y;ti + i°)[h] 

\ / j 1 j — ± 1/7,1 >»3 / v g ^ i — _|_ |^/7,J \ / \ >v\ 


Step 3 


Solve for coefficients 


Coeff's of XT „\ 


Z^-f °( S J z > ^) == E|m| + |n|<2fc ^%-2;m,S S ^ Zm ^ n 


TO 2 

Coeff's of V 2 n ^\q 

2 


N N 

d 2 _ t {s)^F 2 _ 1 {d k _ 2 . mn ;k<^) 


Coeff's of Qjv+i 's 


Qjv-ti (S, Z, z) = E| m | + |n|<W+3 9«±i ;mn (s)2; m ^ 


Off-diagonal Coefficients 


9«+i ;m ,„(0) = -i(l - e^™-)")- 1 - r o ^^ ;min ( S )rf S 


Diagonal Coefficients 


-— L — 

2 ' ^+-;k,k L Jo £t-±-k,k y ' 



Remark Since \fK commutes with complex conjugation, the odd order terms in its Weyl symbol are even 
functions under the involution (s, a, x, £) — ► (s, — a, x, — ^). Although we have restricted to the positive cone 
where a > 0, one gets similar normal forms for a < as long as D 8 is interpreted as \D S \. Since the Harmonic 
oscillators and their powers have even symbols, it must be the case that the QBNF coefficients of the even 
order terms vanish. For instance the coefficient B -2 of must vanish automatically. 

7. Explicit formulae in dimension 2 

To illustrate the algorithm, we carry out the details of the calculation of a 10 , or more importantly the 
normal form coefficients B^.j for k=0, in dimension 2. The result may be summarized as follows: 
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QBNF coefficients for k=0, dim —IThe complete symbol of f (I) in complex coordinates z = y + iij has 
the form B o; 4\z\ 4 + B o . where B o; j are given for both j = 0,4 by weight -2 Fermi-Jacobi polynomials of the 
form: 

Bo,, = \j o + hr\YY\ 2 + b 2T Re(YY) 2 + ct 2 |F| 4 + dr vv \Y\^ + e5 ooT ]ds+ 



0<m,n<3;ro+n=3 IV ' Jl ° 



2 



+ \ Yl C 2;mn Im{[ T u { S )Y m Y n {s)[f T v (t)Y n Y m ]{t)dt]ds}. 

0<m,n<3;m+n=3 " ° " ° 

This corroborates the previous remark that the term B ok vanishes. As will be seen, the corresponding 
density is a total 9 s -derivative. Also, the coefficient Sj of r is the Kronecker symbol, i.e. equals one if j = 
and vanishes if j = 4. 

The expressions for the normal form coefficients in higher dimensions are very similar, and it is only 
for the sake of simplicity that we have stated the result in dimension 2. 

The result for the wave coefficient a 10 is then very similar, modulo the Floquet factors. Indeed, by 
(4.30) or by §5.3 the wave invariant is related to the normal form coefficients by 

a 7 o = T -i{D a ,L)T(a) 

with 



■F„,-i(D a , L) = Lpi(D ai + . . . , D an + L) = L Pl (D ai + i, . . . , D an + 1). 



Hence 



J I 10 2ds = L Pl (D ai + 1 . . . , D an + l)T(a) = L Pl (D ai + !,...,£)„„ + |)H^ =1 (1 - e™* y 1 
In view of (2.33) and the fact that p\ = \Lf a we have 

(7.1) a 70 = J I 1Q ds = \L 2 f (D ai +l,...,D an + |)H7 =1 (1 - e^)" 1 - 

From (7.1) and the formula for the complete symbol of f a one gets the explicit formula for a 7Q as stated in 
the Introduction. 

To prove that the B .j coefficients have the form claimed above, we begin with the expression for f Q 
from (2.33) or from the table in §6: 

(7.2) f (i u . ..,/„) = !/ / v;vlv t \odt. 



o JT"- 



We also have, from (§2, (2.36) and below), that 

(7.3) VI =V + U[Di,Qi}. 

We wish to evalute the coefficients of the complete symbol of f Q in terms of integrals over 7 of Fermi-Jacobi 
data. At first, we will allow the dimension to be arbitrary; when it is time to substitute in metric expressions 
we will restrict to dimension 2. 

Let us consider first the second term on the right side, which simplifies the expression 

[Pi,Qi] + ^ 2 [[i- 1 i? s ,gi],Qi] 

by using [L~ 1 D S7 Qi] = iDi. We recall from §2 or from the table of §6 that the complete symbol of Qi is 
given (in complex vector notation) by 

Qi(s,z,z) = Qi(0,z,z) +L j Vi\ (t,z,z)dt, Qi (0, e M z, e~ ia z) - Qi (0, z, z) = L j Vi\ a (t, z,z)dt. 

Jo J O 
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We note that Vi is independent of D s so that £>i| G = 2? 1. It follows that 

[Z>i (*), Q i (*)] = [©i (a), Q i (0)] + [2?i (a), / © i (*)df] 
so that the second term of (7.3) contributes to / (ii, ■ ■ • , in) the diagonal part of 
(7.4dia 5 ) { [Q i (0, e iQ z, e- fa i),Qi (0, z, z] + 1{[ jf * X>i («)d«, J' X>i (t)dt] . 

Here, the bracket [, ] denotes the commutator of complete symbols in the sense of operator (or complete 
symbol) composition. 

To determine the diagonal part, we write (as usual) 

Q h (0,Z,z)= Yl q h rnn(0)z m z n , Vl( S ,Z,z)= £ di , mn (s)z m z" . 

mn:\m\ + \n\— 3 mn:|m| + |n|=3 

Let us denote the operator composition of two complete Weyl symbols a, b by a#6, that is 

Op w (a) o Op w (b) = Op w (a#b) 
and recall (cf. [Ho III, Theorem 18.5.4]) that this composition has an asymptotic expansion 

i 2 

a#b ~ab + iP x (a, b) + —P 2 (a,b) + . . . 

where Pk(a, b) is the higher order Poisson bracket (or transvectant), a bidifferential operator given in complex 
notation by 

Pk(a,b) = {d z dw - d z d w ) k a(z, z)b(w, w)\ z=w . 
It is well known [loc.cit.] that the commutator is given symbolically by the odd expansion 

a #&-6# a ~ ip 1 ( a ,fe) + -l_P 3 ( a ,6) + ... 

while anticommutators involve only the even transvectants. Since the complete symbols in (7.4 diag) are 
homogeneous polynomials of degree 3, the commutators involve only P\ and P3. One easily computes that 

r\\Z Z ,Z Z ) — lsi ;mnfl „Z Z 

where C\^ mnllv = ^a((m, n), (/j, v)) with a the standard symplectic inner product, and that 

J^3(Z Z ,Z Z ) — ^( m ,n),(p,iy)Z Z 

for certain other coefficients C( m ^^^y A straightforward computation then shows that the diagonal part 
in (7.4 diag) is the sum of the following terms: 

(7.5.1) itf^h (°' e< ^' e ~ ia ~ z )' (°> z < ^ 

+ [C3;303ogi ; 3o(0)gi ; o 3 (0)e l3a +C 3;2 m 
plus 

1 <- L 



(7.5.2) -J J [Di{ S ,z,z),Di{t, Zl z}dsdt 

= ( I [Ci ; 303odl ; 3o( S ) d i;03W+ C l;2112d±^ 
J o J o 

I I [C^Q3Qdl.^(s)d±. m (t) + C^2ll2dl^^ 
J o J o 

We observe that there is no term of order |z| 2 . Using that 

q hmn = - iL (e ia < m -^ lr 1 [ L d hmn ( S )ds 

Jo 

we reduce our problem to the calculation of the di. m?l (t)'s and the d o;m „(t)'s. 
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To evaluate the expressions T>i and V . we conjugate back to the Cs: 

V = ^A* L )C ol i{{A* L )-l, V i = v{A* L )C,_n{Al)-\7.&a) 

with 

, ( LilmY L^ReY 
(7.66) A L ■= 



L-2lmY L-iReY 

Using (2.11) we then compute the Cs in (unsealed) Fermi normal coordinates as follows: 
(7-7.1/2) U=- J2 ^V°M)/ 



3! 

(7-7.1) C = - ]T ^9 oo (s,0)yP- £ d^[g 00 D s + J-^ D s (g oo ji)}\ y=oy 

0:\0\=4 ' 0:\0\=2 

n n 

/3:|/3|=2y=l i,j=l 

+a s 2 + c„r( s ,o) 

where r is the scalar curvature. The last term comes from Ail. 

2 

We then change variables y = Lx and conjugate the symbols. By metaplectic covariance of the Weyl 
calculus, the conjugations change the complete Weyl symbols of the Cs (in the x variables) by the linear 
symplectic transformation Al, i.e by the substitutions 

x -» L~i [(ReY)x + (ImY)Q = \L~^ [Y-z + Y-z] 
' ' £ -> L* [i?eFa; + (/mY)f] = ±£* [Y • z + Yz] 

It is evident that the normal form coefficients are going to be rather lengthy. To give the flavor of the 
full calculations in the simplest setting, we now specialize to the case of surfaces. Later we will briefly extend 
the calculations to all dimensions. 
Dimension 2 

In dimension 2 we have (in scaled Fermi coordinates) 

g 00 (s,y) = l + C 1 T( S )y 2 + C 2 T,(s)y 3 + ... g 11 = 1 J(s, u) = ^gV Q = 1 + C[r( S )y 2 + . . . 

, for some constants Cj,C'j which will change from line to line. Hence the 1/2-density Laplacian in Fermi 
normal coordinates equals 

-A = J-^d^JdsJ- 1 ^ + J-^dyJdyJ- 1 ' 2 

= g°°d 2 + r°d s + d 2 + r^y + g 

and the rescaled Laplacian equals 

-A, = -{Lh)- 2 g$ + 2i(Lh)- 1 g$d s + i^/*) -1 ^ - g$d 2 + T° [h] d s + d 2 y + h~?T\ h] d y + (a) [h] 

Using the Taylor expansion of the metric coefficients one finds that the Cs have the form 

(7.9) £i = CL- 2 T u {s)y\ C a = C x L- 2 y\, v + C 2 L- 1 y 2 T d s + C^r.y 2 - d 2 + C 4 ryd v + C 5 t. 

All terms have weight -2. 

We now switch to the weightless coordinates y = Lx and get: 

(7.9) Ci = CLt„{s)x 3 , C = CiL 2 x a t vv + C 2 L 1 x 2 Td s + C 3 Lt s x 2 - d 2 s + C 4 Txd x + C 5 t. 

Making the linear symplectic substitutions above we first get 

Di(s, z, z) = CL-?t v {s)([Y -z + Y • z]) 3 

hence 

(7-10.1) di. mn (s) = C mn ., 3 L-ir„[Y m ■ Y n ](s) 
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(7.10.2) d hmn (s) = C m n fi L-i r v \Y m ■ Y n ](s)ds 

J o 

(7.10.3) q hmn - -lC mn - 3 (l - e i(m-n)ayl L i J Tu [ Y m ■ Y n }(s)ds 

with m + n = 3 and for certain coefficients C mn;( 3. It is evident that d mn {s) is a Fermi- Jacobi polynomial of 
weight -2 and of Jacobi degree 3, so that the terms (7.5.1-2) are of Jacobi degree 6 as stated in Theorem A. The 
diagonal part has terms of degree |z| 4 and \z\° with coefficients of the form (with m, n = 0, ...3; m + n = 3): 

i fL r-L ~i(n—m)a 

l {[ L TAs)Y m Y«}(s)ds}[J o r y (t)y"y m (t)dt] |1 _ e . (m _ m)a|2 

/L r-L —i(7n—n) 

T v { S )Y n Y m { S )ds][J^ Tu {t)Y m Y n ](t)dt]^-—^-^} 

and 

(7.11.2) j{ f r v {s)Y m Y n {s)[ f T u {t)Y n Y m ](t)dt]ds - f ' T u {s)Y n Y m {s)[ f T v (t)Y m Y n {t)dt}ds}. 

J o Jo Jo Jo 

1 

To calculate the complete symbol of L> s -weight 0, V \ , of the second term Vg we make the same linear 
substitution and eliminate any D s appearing all the way to the right. We also invert the relation 

^Air'D^Al) =D S - \{L- Y dl + Lrx 2 ) 

to get 

^(A^D^A*)- 1 = D S - ^(A^L-'dl + Lrx^niA*)- 1 

and transform the complete symbol of quadratic term by the symplcctic substitution. The result is that 
T> \ (s,z,z) equals 

(7.12.1) dr^lYz + Yz} 4 

(7.12.2) +C 2 L- 2 t[Yz + Yz] 2 #{[Yz + Yz} 2 + tL- 2 [Yz + Yz} 2 )) 

(7.12.3) +C 3 t s [Yz + Yz} 2 

(7.12.4.1) -2L- 2 d s ([Y z + Yz} 2 - L- 2 t[Yz + Yz] 2 ) 

(7.12.4.2) +L- 2 {[Y z + Yz} 2 - L~ 2 t[Y ■ z + Yz} 2 }#{[Yz + Yz} 2 - L~ 2 t[Yz + Yz} 2 } 

(7.12.5) +C i L- 1 T(Yz + Yz)#{Yz + Yz) + C 5 t. 

Our concern is with the diagonal part of the complete symbol, that is, with the terms involving 
|z| 4 , \z\ 2 , \z\°, and more precisely with their integrals over 7. To begin with, we observe that the diago- 
nal part of term (7.12.1) is purely of degree \z\ 4 and its average over 7 equals 

1 ' L 



(7.13.1) (Const.)\z\' 1 ■ - j r^Y^ds. 

The diagonal part of term (17.12.2) contributes only the P and Pi terms, of degrees |z| 4 and \z\° respectively, 
whose averages over 7 have the form 

(7.13.2) {\z\ 4 /or/\z\°) ■ \ f ' T\ar\Y\ 4 + bRe{YY) 2 + c\YY\ 2 ]ds 

L J 

with constants a, b, c, d which can differ between the two degrees. The missing Pi-term vanishes: it is a 
multiple of the Poisson bracket 

PidYz + Yzf^iYz + Yz} 2 )) 
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which simplifies to a term of the form 

t[Y 2 Y 2 - Y 2 Y 2 } = t(YY - YY)(YY + YY) = Ct-^-\Y\ 2 

by the symplcctic normalization of the Jacobi eigenfield. However the integral over 7 of this term vanishes, 
that is 

(7.13.3) f t s \Y\ 2 =0. 

J o 

This can be seen from the Jacobi equation, which implies: 

[Y(Y')" + t s \Y\ 2 + tY'Y] = 0; 

integrating over 7 and integrating the first term by parts twice kills the outer terms and hence the inner 
one. In a similar way, the diagonal part of (7.12.3) is of pure degree \z\ 2 with coefficient r s |F| 2 , so it 
makes no contribution either. Nor does the term (7.12.4.1), which is manifestly a total derivative and hence 
automatically has zero integral. The term (7.12.4.2) is a composition square and hence contributes only a 
product P c -term of degree \z\ 4 and a P2-term of degree \z\°, namely (for j=0,2) the diagonal part of 

Pj(z 2 Y '' + 2\z\ 2 \Y\ 2 + z 2 Y 2 + t(z 2 Y 2 + 2\z\ 2 \Y\ 2 + z 2 Y 2 , z 2 f 2 + 2\z\ 2 \Y\ 2 + z 2 Y 2 + t(z 2 Y 2 + 2\z\ 2 \Y\ 2 + z 2 Y 2 ) 
whose average over 7 has the form 

(7.13.4.2) {\z\ 4 /or/\z\°) ■ \ / L [a|F| 4 + bTRe{YY 2 ) + cr\YY\ 2 + d T 2 \Y\ A ]ds 

L J 

where again the coefficients may vary between the two degrees. Finally, we the first term of (7.12.5) obviously 
has no diagonal part while obviously the second term contributes 

(7.13.5) c\ I rds. 

L Jo 

This completes the analysis of the QBNF coefficients B 4, B 2, B a o- 

8. Appendix : The classical Birkhoff normal form 

The method of §2 for putting \fK into quantum Birkhoff normal form began, essentially, by putting 
the linearization of >/~A at 7 into normal form by a linear symplectic transformation, and then proceeded by 
induction on the jet filtration at 7. The purpose of this appendix is to describe, rather briefly and sketchily, 
how to put the principal symbol of y/A into Birkhoff normal form by an analogous method. (We have not 
found this particular algorithm in the literature, but it is quite likely that it, or a much better algorithm, is 
well-known). We hope that the classical algorithm will help clarify the procedure in the quantum case. 

In the usual Fermi symplcctic normal coordinates (s, a, y, r/), we may write the principal symbol of vA 
in the form 

n 

H(s,a,y,r]) := (g 00 (s, y)a 2 + ^ g v \s, y)ViVi) h 

ij=l 

n 

tj=l 

Taking the Taylor expansion at y = rj = we get 

H(s, a, y, n) - (7(1 + KiMViVi + E ^ + • • • ) 

ij—i i—i 

from which we extract the linearized symbol 

1 " " n 2 

(A.l) h{s, a iy ,ii)=<j + -J2 K ij {s)yiVi° + ~f- 

ij=i 1=1 a 
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We make the symplectic change of variables 

1 yvj fj 

(A.2) 4> : (s,a,y,r}) -> (s' ,a' ,y' ',?/) := (s + - — ,^y,—=) 

I a ^Ja 

which transforms h into 

^ n n 

(A3) h{s, a,y,7]) = a + - ^ K l0 {s)y iyj + ^ r}\ 

ij — l i—1 

and H into 

H{s,a,y,rj) = h + h [3] + ... 
with . . . denoting terms of order 3 and higher in (y, rf). (Such terms are of order 3/2 with respect to the order 
in the isotropic calculus, while h is of order 1, which is the rationale for calling h the linearized symbol). 
The first step in putting H into Birkhoff normal form is to put h into Birkhoff normal form 



1 ™ 



2 U-. 



by means of a symplectic map. Equivalently, we wish to convert the Hamilton equations 

-fs = 1 

as 

-fo = 

as 



£v = -K{s)y 



as 



into the linear equations 

Ts a = a P 
TsP = -°9- 

We first do this in just the (y, rj, q,p) variables, with a symplectic map of the form (y, 77) — > (q,p) = B(s){y, r)). 
The condition on B is that 

(AA) BB- 1 + BKB- 1 = aj 

where K is the block anti-diagonal matrix with blocks I and —K, where J is the usual block anti-diagonal 
matrix with blocks ±7 and where a J is the block anti-diagonal matrix with coefficients ±aj replacing ±1 
in J. The equation (A. 4) is equivalent to 

(A5) -^-B- 1 +K{s)B- 1 aB- 1 J 

as 

which has the solution 

(A6) B(s)=r a (s)a- 1 . 

To make the map symplectic in all the (s, cr, y, rj) variables, we observe that the map 

1 " 

i>\ ■ {s,v,y,v) ^{s,a- 2^«i(j/i +Vi),r a {s)(y,r))) 

is symplectic with respect to ds A da + dy A drj and satisfies 

1 ™ 

^1 ( CT + 9 I>*(^ 2 = CT - 
1 i=i 

Also, the map 

■02 : (s, a, y, 77) -> (s, a + f(s, y, ij),a~ 1 (y, r?)) 
with f(s,y,r]) :~ — i X^=i ( s )yiUj + S"=i ?7? i s symplectic with respect to ds Ada + dyA dr\ and satisfies 

i)\(a) = a + f. 

It follows that xi := 02 • 01 pulls back /i to its Birkhoff normal form. 
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Let us now write 



Hi ~ x * 1 (H) = h + hf ] +hW + . 



with hf^ (s, a, y, rf) of order 1 and vanishing to order k at (y, jf) = (0, 0). Following the algorithm for putting 
a Hamiltonian with non-degenerate minimum at into Birkhoff normal form [AM, 5.6.8, p. 500] we seek a 
symplectic map of the form 

X2 = expada? F$ 
with F3 = Fz(s,y,rj) vanishing to order 3 and with 

(A.7) {h, [x^iFl = 0. 

Here, {, } denotes the Poisson bracket with respect to ds A da + dy A dr\ , [F]- k denotes the terms of vanishing 
order < k and adF ■ G := {F, G}. The equation (A.7) is equivalent to 



(A.8) 



hf ] + {a?F 3 , h} e ker(adh) n V 3 



where V k denotes the space of homogeneous polynomials of degree k in (y, rj) with smooth coefficients in s 
and with an overall factor of a 2 _1 . In terms of the complex coordinates Zj = yj + ir]j on R 2n we may write 

b1+|fe|=3 



with j = (j 1 ,...,j n ),k= (fci, 
vector field of h as 



E a As)z j z k 

|j| + |fe|=3 

,k n ), and we may write the Lie derivative C with respect to the Hamilton 

d . d _ d . 



dzi 



&Zi 



The monomials in keradC n T 3 * are of the form z 3 z with (a,j — k) = which implies j 
assumptions on a. No such terms occur for odd k, so equation (A.8) thus becomes 



k with our 



(A.9) 



bl+|fc|=3 



k(s) + a ■ (j - k)c jk (s) = -a jk (s). 



Expanding the periodic (or more generally almost periodic) functions Cjk and ajk in Fourier series 



c 3k (s) = e lms c jk (m) a jk (s) = ^ e vms a ]k (m) 



we can solve (A.9) with 
(A10) 



Cjfc(m) 



a jk {m) 



im + (a, j — k) 



Writing H2 — X2 ' X\H , we arrive at the analogous problem for the fourth order terms. As in the 
quantum case, the even steps behave a little differently from the odd ones since now there can be terms 
in f/J, 4 ' with \j\ = \k\ and hence which lie in keradC n "P 4 . Since these terms already commute with L it 
suffices to solve the analogue of (A9) with only the coefficients a,j k (s) coming from H% ] - H^ ]o . 

This puts the terms up to fourth order in normal form, and the process continues inductively to define 
symplectic maps \n ■ Xn- 1 Xi which pulls back H to a normal form up to degree N. Since \n — ImodOjss 
the infinite product defines a smooth symplectic map which pulls back H to a normal form modulo Ooo, 
that is, to its Birkhoff normal form. 
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9. Index of Notation 



In the following, tl will denote the translation operator tl/(s, y) = f(s + L, y) on functions on K. x R n . 
Operators A then transform under tl by t^At* l . 
NI.l: Model objects 



Model object 




Functions 


tlJ = f 


Vy L/V_^i. U LUl il) 


Tr At* — A 


Maximal abclian algebra 


A =< 1Z, h J„ > 


Distinguished clement 


TI := ±(LD S + H a ) 


Harmonic Oscillator 


H a := | Efc=i a kh 


Gaussians/Hcrmites 


lo(y) = liM = lq = C q A*^ . . . A* n *<* 


,4-Eigcnfunctions 


<t>l q {s,y) ■= efe(s) ®7g(y) 


"^-Eigenvalues 


K4>° kq {s,y) = r kq <P° kq {s,y), r kq = ±{2-Kk + E"=i(?j + l) a i) 


Scaled Laplacian 


TZ h := /i(a)A h /i(5)* ~ E™ =0 h^+^n^ 


Intertwiner to SC normal form 


W h := n^W^^a), := exp{ih?Q,) 


S.C. Normal form 


W^n h W h \ ~ /i" 2 + E'^o Wfc(Ji, ■ ■ ■ , fn) 



NI.2: Twisted model objects 



Twisted model object 




Functions 


r L f = n(r a {L))f 


Operators 


t l At* l = »(r a (L))A»(r a (L))* 


Maximal abelian algebra 


A=< D s ,h,...,I n > 


Distinguished clement 


D s 


Harmonic Oscillator 


H a := J Efc=l a k l k 


Gaussians /Hcrmites 


lo{y) = 7i/(9) = c-*i^ 79 - O^! 91 . . . a;*- 


„4-Eigenfunctions 




D s -Eigenvalues 


D s e"' k " s <g> 7,(y) = r kq e lrkl > s <g> 7 9 (j/) 


Scaled Laplacian 


% := Ma)A fc Ma)* ~ E™= /* ( " 2+ ^2-^ 


Intertwiner to SC normal form 


Wfc ~ ^=iW h ^(r a ), W h t := exp(ih%Q±) 


S.C. Normal form 


w* h v h w h \ ~ ft- 2 + e^o fc fc /fc(A, ■ ■ • , 4) 



NI.3: Adapted model objects 



Adapted model object 




Functions 


tlS = KT).f 


Operators 


r L Arl = fi(T)A(j,(T)* 


Maximal abelian algebra 


A =< C, AiAJ, . . . , A„A* > 


Distinguished element 


C = D S - i(E -=i ^ + E«=i 


Harmonic Oscillator 


ff« := 5 ELi "feAfeA* 


Gaussians /Hcrmites 


MOloM := t/ (s,u) = (deiF^))"^^! < r(«)u,u >, T( S ) := ^F" 1 ; 
/i(a" 1 )7 9 :=t^ = C g Af 1 ...A«-C^. 


*4-Eigenfunctions 


M5.XO = fo, :=e"^£/ g ( S ,u) 


Scaled Laplacian 





NI.4 Intertwining operators 

1. /z(r a ) := J gl fi(r a (s))ds 

2.a s :=( i m ^ S ): l m 5 S l )=A 



3. ^l(s) := 



i?eT(s)* i?eT(s)* 

LilmY(s) L~i ReY(s) 
L-ilmY(s) L-^ReY(s) 
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Intertwining operator A 


t l At£ 


fj,(a) : Adapted Model -» Twisted Model 


tl^{o)t* l = 


fi(r a )n(a) 


fj,(r a ) : Model -> Twisted Model 


T L ^{r a )T* L -- 


= n{r a {L))n{r a ) 


With a s := a s r a (s), fj,(a): Adpated Model — > Model 


rMa)r* L = 


fj,(a) 
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